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Abstract. The covariant density functional theory with a few number of parameters 
has been widely used to describe the ground-state and excited-state properties for 
the nuclei all over the nuclear chart. In order to describe exotic properties of 
unstable nuclei, the contribution of the continuum and its coupling with bound states 
should be treated properly. In this Topical Review, the development of the covariant 
density functional theory in continuum will be introduced, including the relativistic 
continuum Hartree-Bogoliubov theory, the relativistic Hartree-Fock-Bogoliubov theory 
in continuum, and the deformed relativistic Hartree-Bogoliubov theory in continuum. 
Then the descriptions and predictions of the neutron halo phenomena in both spherical 
and deformed nuclei will be reviewed. The diffuseness of the nuclear potentials, nuclear 
shapes and density distributions, and the impact of the pairing correlations on nuclear 
size will be discussed. 
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Summary and perspectives 


1. Introduction 


The building blocks of nucleus are protons and neutrons. An arbitrary combination of 


protons and neutrons doesn’t necessarily make a nucleus. Among thousands of nuclides 


expected to exist, only 288 of them are stable or at least with half-lives longer than the 


expected life of the Solar System, which form the valley of stability in the nuclear chart. 


These nuclides are sandwiched by about 3,000 nuclides produced in laboratories [1]. 
The nuclear density functional theory (DFT) has predicted that about 7,000 nuclides are 
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bound with respect to neutron or proton emission |2, 3, 4]. If the continuum contribution 
is included, as demonstrated in Ref. [5] from O to Ti isotopes, this number may exceed 
10,000. In addition, there are also nuclides (resonances) which, though beyond the drip 
lines, can have lifetimes long enough to be measured and investigated. 

The nuclides lying above the valley of stability in the nuclear chart are proton-rich 
and the proton drip line has been identified experimentally up to protactinium (Z = 91). 
However, for those lying below the valley of stability, which are neutron-rich, the border 
of the nuclear territory is known only up to oxygen (Z = 8). Experimental exploration of 
very neutron-rich nuclei is extremely challenging because of very low production rates 
in studies involving the fragmentation of stable nuclei, as well as the separation and 
identification of the products. 

In nuclear physics, the study of the properties of exotic nuclei — nuclei with extreme 
numbers of proton or neutron — is one of the top priorities as it can lead to new insights 
on the origin of chemical elements in stars and star explosions. Although the radioactive 
ion beams (RIB) have extended our knowledge of nuclear physics from stable nuclei to 
exotic ones far away from the valley of stability, it is still a dream to reach the neutron 
drip line up to mass number A ~ 100 with the new generation RIB facilities developed 
around the world, including the RIKEN Radioactive Ion Beam Factory (RIBF) in 
Japan |6], Facility for Antiproton and Ion Research (FAIR) in Germany [7], Second 
Generation System On-Line Production of Radioactive Ions (SPIRAL2) at GANIL in 
France [8], Facility for Rare Isotope Beams (FRIB) in USA [9], Rare Isotope Science 
Project (RISP) in Korea [10], Cooler Storage Ring (CSR) at Heavy Ion Research Facility 
in Lanzhou (HIRFL) in China [11, 12], etc. 

The development of RIB facilities around the world significantly stimulates the 
study of exotic nuclei far from the valley of stability [13, 14, 15, 16, 17, 18, 19, 20, 21, 
22, 23]. New and exotic phenomena were observed in nuclei close to drip lines such 
as neutron or proton halos [24], changes of nuclear magic numbers [25], the island of 
inversion [26], pygmy resonances [27], etc. More exotic nuclear phenomena have also 
been predicted, e.g., giant halos [28, 29], shape decoupling between core and halo [30, 31], 
etc. 

In halo nuclei, the extremely weak binding of valence nucleons leads to many new 
features, such as the coupling between bound states and the continuum due to pairing 
correlations and the very extended spatial density distributions. Therefore, one must 
consider properly the asymptotic behavior of nuclear densities at a large distance from 
the center and treat in a self-consistent way the discrete bound states, the continuum, 
and the coupling between them in order to give a proper theoretical description of the 
halo phenomenon [20, 32]. There are mainly two key points to achieve a self-consistent 
description of halo nuclei with density functional theories: the pairing correlations and 
the suitable basis. 

In the mean field descriptions of nuclear halos, nuclear pairing plays an important 
role [33, 34, 35, 28, 29]. The conventional BCS method for the pairing is not appropriate 
for exotic nuclei. To include properly the contribution of continuum states, the 
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Bogoliubov transformation has been justified to be very useful [36, 37]. By solving 
the nonrelativistic Hartree-Fock-Bogoliubov (HFB) [36, 37] or the relativistic Hartree- 
Bogoliubov (RHB) [34, 38, 39] equations in coordinate (r) space, one can fully take into 
account the coupling to the continuum. As an alternative, one may first locate resonant 
states in the continuum and compute the resonance parameters with various methods 
[40, 41, 42, 43, 44, 45, 46], then by using the resonant-BCS (rBCS) approach [47], include 
the contribution of these resonant states and study halo phenomena [48, 49, 50]. The 
Green’s function method [51], which can take into account properly the asymptotic 
behavior of continuum states, has also been adopted for solving the HFB [52, 53, 54], 
the self-consistent Skyrme-HFB [55, 56, 57] or RMF [58] equations in the coordinate 
space. 

The solution of the coupled differential equations of HFB and RHB theories is 
relatively easier for spherical systems with local potentials, where one-dimensional 
Numerov or Runge-Kutta methods can be applied. This is also true even for nonlocal 
problems where the finite element method (FEM) can be used [39]. A more widely used 
method for solving such equations is to expand the single particle wave functions in a 
basis and this is more convenient particularly for deformed systems. In the past decades, 
the harmonic oscillator (HO) basis has been used with a great success for stable nuclei 
in both nonrelativistic and relativistic mean field approximations such as the Skyrme 
Hartree-Fock, Hartree-Fock-Bogoliubov, relativistic Hartree, and relativistic Hartree- 
Bogoliubov theories for deformed or nonlocal systems [59, 60, 61]. However, due to the 
incorrect asymptotic behavior of the HO wave functions, the expansion in a conventional 
HO basis is not proper for the description of exotic nuclei [62]. Considerable efforts 
have been made to develop mean field models in a basis with an improved asymptotic 
behavior at large distances, e.g., a transformed HO basis [63, 64] or a Woods-Saxon 
(WS) basis [65]. 

The DFT with a small number of parameters is widely used to study the ground- 
state and excited-state properties of the nuclei all over the nuclear chart. In particular, 
the covariant density functional theory (CDFT), incorporating the Lorentz symmetry 
in a self-consistent way, has received wide attention [66, 67, 68, 20, 69, 70]. 

There exist a number of attractive features to describe atomic nuclei in the 
relativistic framework. The most obvious one is the natural inclusion of the nucleon 
spin degree of freedom, resulting in the nuclear spin-orbit potential that emerges 
automatically with the empirical strength in a covariant way. Moreover, it provides 
a new saturation mechanism for nuclear matter [71], reproduces well the measured 
isotopic shifts in the Pb region [72], reveals more naturally the origin of the pseudospin 
symmetry [73, 74] as a relativistic symmetry [75, 76, 77, 78, 79, 80, 81, 82], and predicts 
the spin symmetry in the anti-nucleon spectrum [83, 84]. It can also include the nuclear 
magnetism [85], i.e., a consistent description of currents and time-odd fields, which plays 
an important role in the nuclear magnetic moments [86, 87, 88, 89] and nuclear rotations 
(90, 91, 92, 93, 70], etc. 

The CDFT is a reliable and useful model for nuclear structure study in the whole 
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nuclear chart, and therefore, it is natural and necessary to describe the halo nuclei based 
on the CDFT. 

Since the discovery of the neutron halo, i.e., the anomalously large nuclear radius in 
“Ti [24] which deviates from the conventional A'/? law on the mass dependence of the 
nuclear radius, considerable efforts have been undertaken to understand this interesting 
phenomenon in exotic nuclei close to the neutron drip line. 

The relativistic mean field models have been used to study nuclear halo phenomema 
since early 1990s [94, 95, 96, 97, 98, 99, 100]. By incorporating the Bogoliubov 
transformation in the relativistic Hartree theory, the relativistic continuum Hatree- 
Bogoliubov (RCHB) theory has been developed and extensively used to study the 
halo phenomeon in spherical nuclei [34, 38, 28, 101, 39]. The RCHB theory was also 
extended to including the hyperon and used to study neutron halo in A-hypernuclei 
[102, 103, 104, 105, 106]. The relativistic Hartree-Fock-Bogoliubov theory in continuum 
(RHFBc) theory has been developed and used to study the influence of exchange 
terms on the formation of neutron halos [107]. Over the past years, lots of efforts 
have been made to develop a deformed relativistic Hartree-Bogoliubov theory in 
continuum (DRHBc) [108]. Halo phenomena in deformed Ne and Mg isotopes have 
been investigated with the DRHBc theory [30, 31, 109]. 

In this Topical Review, recent progresses on the application of the CDFT for the 
neutron halo phenomena in spherical and deformed nuclei will be presented. As an 
extensive review is available [20], the RCHB theory and its application will be briefly 
mentioned here for completeness. In Section 2, we give the formalism for CDFT in 
continuum, including the general framework of CDFT, the formalism for RHB theory, 
and the RCHB, RHFBc and DRHBc theories. In Section 3, the applications of the 
CDFT for halos and giant halos in spherical nuclei will be reviewd. We will discuss the 
RCHB description of the neutron halo in !'Li, the prediction of giant halos in Zr and Ca 
isotopes, the continuum contribution and the extension of nuclear landscape, halos from 
relativistic and non-relativistic approaches, and halos with RHFBc theory. The densities 
and potentials in exotic nuclei become more diffuse and the surface diffuseness has 
important influences on spin-orbit splittings, single particle level structure, and nuclear 
masses. Recent progresses on these topics will be presented in Section 4. In Section 5, 
we will discuss the influence of pairing correlations on the nuclear size. The necessity to 
carry out a self-consistent calculation will be shown in order to understand the impact 
of the pairing correlations on the nuclear size. In Section 6, the study of neutron halo in 
deformed nuclei will be reviewed, including the neutron separation energies, quadrupole 
shapes, density distributions, and radii. The generic conditions for the occurrence of 
halos in deformed nuclei and the predicted shape decoupling between the core and the 
halo will be discussed. Finally we give a summary and discuss perspectives in Section 7. 


201707.00914v1 


chinaXiv 


CONTENTS 6 


2. Covariant density functional theory in continuum 


2.1. General framework of CDFT 


The covariant density functionals represent a number of relativistic mean field (RMF) 
and relativistic Hartree-Fock (RHF) models based on the framework of quantum 
hadrodynamics (QHD) [66]. The advantages in using covariant functionals include the 
natural inclusion of the nucleon spin degree of freedom, the consistent treatment of 
isoscalar Lorentz scalar and vector self-energies which provides a unique parametrization 
of time-odd components of the nuclear mean field, a natural explanation of the empirical 
pseudospin symmetry, and a new saturation mechanism of nuclear matter with a 
distinction between scalar and four-vector nucleon self-energies [67, 110, 111, 68, 20, 
112, 82]. 

The conventional CDFTs are based on the finite-range meson-exchange 
representation, in which a nucleus is described as a self-bound system of Dirac nucleons 
interacting with each other via the exchange of mesons through an effective Lagrangian 
density, which is associated with the nucleon (Y), the o-, w-, p-, and 7-meson fields, and 
the photon field (A), 


L = Ly + Lo + Lou + Lo + Lr + Lat Lr, (1) 
containing the Lagrangian densities for the nucleon Ly, 

Ly =} (iy, — M) Y, (2) 
for the o-, w-, p-, and 7-mesons, and the photon A, 


1 1 
Lo = + z oOo = Moos 


1 1 
£ = = g Ow + sone, 
1 > w, l 2 > 
L= — giw R" + ZPP * Pus 
1 
Lr = + aa - OLT — mk - T, 
I py 
La = = i P ii 


and for the interactions between nucleon and meson fields, 


7 sos Jo mas s 
Lr =Ņ |90 — Ju Y Wy = Ora * Pu T T Oy Pu eT 


tn > oS f= T3 
=m WOT -T — eg" J Ay| Y. (3) 


In above Lagrangian densities, M denotes the mass of nucleon, and m, (go), Mw 
(gu); Mp (gp, fp), and My (fr) are respectively the masses (coupling constants) for o-, w-, 
p-, and z-mesons. The p-tensor coupling term, which is practically negligible at Hartree 
approximation, can significantly improves the descriptions of nuclear shell structures at 
Hartree-Fock approximation [113]. 
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The field tensors of the vector mesons Q and R and the electromagnetic field F are 

defined as, 
Q = Fw — ow”, 
Re = apr — arm, (4) 
F” = oh A” — OV A". 

The isoscalar-scalar g-meson provides the long-range attractive part of the nuclear 
interaction whereas the short-range repulsion is governed by the isoscalar-vector w- 
meson. The photon field A“(2) accounts for the Coulomb interaction and the isospin 
dependence of the nuclear force is described by the isovector-vector p-meson. The m- 
meson field and the rho-tensor coupling contribute only if the Fock terms are included. 

In the following, the vectors in the isospin space are denoted by arrows and the 
space vectors by bold type. Greek indices u and v run over the Minkowski indices 0, 1, 
2,3 ort, x, y, z, while Roman indices 7, j, etc. denote the spatial components. 

The variation of the Lagrangian L = f d’xL£(x) with respect to nucleon (Y) and 
meson fields (¢;) leads to the Dirac equation of nucleon, 


(=07"0, + Md) pla) =0, (5) 


with the self-energy = to be self-consistently determined, and the Klein-Gordon 
equations for the meson fields (¢ = o, w,,, Pu and 7), 


("O, +m) $ = So, (6) 


with the source terms 


—gobh, o=o, 
5, - IVW, f $ = Wy, E 
IPT +P r owh, $= Pu, 
FP dys w TY, =F, 
and the Proca equation for the photon field (A), 
O Fyn = ep aa. (8) 


Under the no-sea and mean field approximations, the Dirac Hartree-Fock equation 
can be derived as 


f dr'h(r, r')b(r!) = eb (1), (9) 


where € is the single particle energy (including the rest mass) and the single particle 
Dirac Hamiltonian h(r,r’) contains the kinetic energy h™, the direct local potential 
hP, and the exchange non-local potential hE, 


h®”(p, 9!) = [ov p+ BM] alr — (10) 
h?” (r, r’) = [Er(r)ys + Do(r) + BUs(r)] 6(r — r’), (11) 


i "E Yelrir) Yp(r,r’) 
h®(r,r’) = Cn Lie l (12) 
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In the above expressions, the local scalar Xs, vector “io, and tensor Xr self-energies 
contain the contributions from the direct (Hartree) terms and the rearrangement terms, 
while the nonlocal self-energies Xa, Xr, Yo, and Yr come from the exchange (Fock) 
terms. More details can be found in Ref. [114]. 

If one sticks to the Hartree level (RMF), for a nucleus with time reversal symmetry, 
the Dirac equation (5) is reduced as 


h(n) = evs(r), (13) 


where 


A 


h=a-p+6[M +%Xs(r)] + Xo(r), (14) 
is the single particle Hamiltonian. 

More recently, the RMF framework has been reinterpreted by the relativistic Kohn- 
Sham density functional theory, and the functionals have been developed based on the 
zero-range point-coupling interaction [115, 116, 117, 118]. In fact, by solving formally 
the Klein-Gordon equations for mesons and taking the large meson mass limit, one can 
derive point-coupling functionals in which the meson in each channel (scalar-isoscalar, 
vector-isoscalar, scalar-isovector, and vector-isovector) is replaced by the corresponding 
local four-point (contact) interaction between nucleons. 

The basic building blocks of CDFT with point-coupling vertices are 


(POTY), O = {1,7}, re {1, Yu Y5, Ve tas Tuv t, (15) 


where w is the Dirac spinor field of nucleon, 7 is the isospin Pauli matrix, and T 
generally denotes the 4 x 4 Dirac matrices. There are ten of such building blocks 
characterized by their transformation characteristics in isospin and Minkowski space. 
A general effective Lagrangian can be written as a power series in YOT% and their 
derivatives, with higher-order terms representing in-medium many-body correlations. 
For details, see Refs. [112, 115, 116, 117, 118]. 

The point-coupling model has attracted more and more attention owing to the 
following advantages. First, it avoids the possible physical constrains introduced by the 
explicit usage of the Klein-Gordon equation to describe the mean meson fields, especially 
the fictitious ø meson. Second, it is possible to study the role of naturalness [119, 120] 
in effective theories for nuclear structure related problems. Third, it is relatively easy 
to include the Fock terms [121], and provides more opportunities to investigate its 
relationship to the nonrelativistic approaches [122]. 


2.2. Relativistic Hartree-Fock-Bogoliubou theory 


Pairing correlations are crucial in the description of open-shell nuclei. For exotic nuclei, 
the conventional BCS approach turns out to be only a poor approximation [37, 20]. 
Starting from the RMF Lagrangian density, a relativistic theory of pairing correlations 
in nuclei has been developed by Kucharek and Ring [123, 67], 


J He.r’) ( vu(r’) p ( pur) | (16) 
by (r’) vv (r) 
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with 
h(r,r’) — A}0(r — r’ A(r,r’) 
H(r, r’) = ) , (17) 
A(r,r’) —h(r, r’) + Aor = r’) 

where the chemical potential A is introduced to preserve the particle number on the 
average. 

The single particle Hamiltonian h(r,r’) is derived from Eq. (5) with the 
retardation effects neglected, e.g., the Hamiltonian given in Eq. (13) under the Hartree 
approximation. The pairing potential can be written as, 


Aal n= ae kg(r, r’), (18) 


where the pairing tensor «K is 
Kalt, r’) = Yva (r)" Yua (T). (19) 


The pairing force is either taken as a density-dependent two-body force in a zero range 


limit, 
A 1 1 1 p(r) 
V(r,r’) = V(r — r')—-(l-a@- 0’) ( 1-—}], (20) 
4 Po 
with an adjusted strength Vo, or as the pairing part of the Gogny force [124], 
V(r, r’) = So el 1/4)" (W, + BiP? — HP — M,P7P’), (21) 
i=1,2 


with the parameters pi, Wi, Bi, Hi and M;. It has been shown that after proper 
renormalization, the above density-dependent force of zero range and the finite range 
Gogny force produce essentially the same results [35]. 


2.8. Relativistic continuum Hartree-Bogoliubov theory 


If the Fock terms are neglected, as it is usually done in the covariant density functional 
theory, the Dirac Hartree-Bogoliubov equation for the nucleons still takes the same form 
as that given in Eq. (16). For spherical nuclei, the quasiparticle wave functions read, 


= ean | a 
am Ern) “i 


The angular part of the RHB equation can be integrated out and the radial RHB 
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equations are derived as 
dGy(r) œ% 


UNI + EGulr) — (E+ = Zolr) + Bs(r))Fulr) +r fr dr Ap(r,r)Fy(r) = 
FH) _ ER) + (B+A— Solr) —Bslr))Gutr) +r f rarer Gv"). = 
ge) + =Gy(r) + (B= A+ Zo(r) — Bsr) Fv(r) +r f r'dr'Ap(r,r')Fy(r') = 
EED L By (0) — (E -A+ Eole) + Bs()Gv(r) +r | rdr Aatrau) = 


The pairing potentials are ae as 


Ag(r,r’) m "rg X GEGE) + GENGE), (25) 


Eg>0 


o ie Noe SEG) Fv (r) + FE) FE). (26) 
E,>0 
If the zero-range pairing laten is used, the above coupled integro-differential equations 


are reduced to differential ones, which can be directly solved in coordinate space [125]. 


2.4. Relativistic Hartree-Fock-Bogoliubov theory in continuum 


For spherical systems, the solution of the RHFB equations, i.e., the Dirac spinor Yy, 
and wy,, can be written similarly to Eqs. (22) and (23), 


1 IGu) YVE ml) 
pu, (r) = - 21 
a am al a 
1 Ov) ng (*) 
Yy (r) = — 28 
Gi a m 


The RHFB equations (16) are then reduced to the coupled integro-differential equations, 


= siar Gu, (r) — (Ea + A- E) Fy, (r) 


dr 
+ Xy,(r) + r | rar’ Bol ely le) =0, (29) 
ds Ka 
[E-E - Dr) Posl) + (Ea +A- EB) Gule) 


Ta ; ey / rdr'A,(r,r')Gy, (1) =0, 
= +24 Er Gy, (r) + (Ea —A+Z_) Fy, (r) 

ore | rdr'A,(r,r")Fu,(r) = 0, 
[E-E Br] Fal) (Ea -A+ EN Gul 


—Yy,(r)+r | r'dr'Ag(r,r')Gu,(r’) = 0, 
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where E, are the quasi-particle energies (without the rest mass), and the local self- 
energies X, and %_ are 


by = Xo + Us, (30) 
O_ =D) — Es — 2M. (31) 


More details can be found in Ref. [114]. 
The pairing potentials A,(r,7r’) in Eq. (29) can be expressed as 
Avr ==) VP Or yale); (32) 
b 


where the pairing tensor «(r,7r’) reads 
1. i” 
nalr r’) = 5 [Gu (r)Gu (r) + Fal) Fv (r+ [Gv (r)Gu, (r) + Fv (r) Ful’), (33) 


where j? = 2j +1. Details of the pairing interaction matrix element V4? can be found 
in Ref. [39]. 

In contrast to the RHB approach with 6-force in the pairing channel where the 
radial equations (29) become differential equations, in RHFB theory the radial equations 
are fully integro-differential. The integral terms may arise from the Fock terms or the 
pairing channel if the finite-range pairing force is used. In coordinate space, it is difficult 
to solve such equations by the localization procedure adopted in Refs. [126, 127] when 
solving the relativistic Hartree-Fock equations. It is more convenient to solve them by an 
expansion of the Dirac-Bogoliubov spinors in an appropriate basis. This has been done 
by using the Dirac Woods-Saxon (DWS) basis [65]. This basis has been constructed for 
the investigation of weakly-bound nuclei. The set of DWS basis functions 


{les 9a(7,7)] ; € Z OF, (34) 


are eigenfunctions (with eigenvalues ¢,) of a Dirac equation with Woods-Saxon-like 
potentials for Mo(r)+Ns(r). They are determined by the shooting method in coordinate 
space within a spherical box of size Rmax [128]. 

The U and V components of the Dirac-Bogoliubov spinors (27) can be expanded 
as, 


Nr Np 
bu = >) Upp + > Uiga, (35) 
p=1 d=1 


Np Np 
py = X "Ved + X Vaga, 
p=1 d=1 


where Np and Np respectively correspond to the numbers of positive (ep > 0) and 
negative (£q < 0) energy states in the DWS basis. Obviously, because of spherical 
symmetry the quantum number « is preserved, i.e., the RHFB equations can be solved 
for each value of « and the sums in the expansion (35) run only over states with the 
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same «K. For a fixed value of x we have the radial basis spinors 


(36) 
(37) 


where the sub-indices p and d correspond to the number of nodes of the basis functions 
Gp for positive energy and Fy for negative energy. 

In the DWS basis (35) the radial RHFB equations (29) are transformed to a matrix 
eigenvalue problem, 


TEn a. ? 


where H and A are (Nr + Np) x (Nr + Np)-dimensional matrices, U and V are the 
column vectors with Np + Np elements. From the expressions of the single particle 
Hamiltonian A and pairing potential A given in the previous part we obtain the matrix 
elements of H and A as 


Hen = f drGn(r) (-4 $ £) Fy(r) + f dr F, (r) (= ie £) Gw (rX39) 
7, = f dr [Galt Gar (r)E.(r) ARORO 
f dr|Ga(r)Fw(r) + Ga(r)Fw(r)] U(r), (40) 
"E , Yo Yr Gy 
Huw = Jefe (Gn Pa), ( Xo Xp l ( Fy E an 
App = fo | arate, r’) (GaGn(r’) + Fi(r) Fa (r’)] , (42) 


where n and n’ run over the radial quantum numbers of the DWS basis states in Eq. (37) 
with both positive energies (n, n’ = p) and negative energies (n, n’ = d). 


2.5. Deformed relativistic Hartree-Bogoliubov theory in continuum 


For deformed nuclei the solution of HFB or RHB equations in r space is numerically 
a very demanding task. Considerable efforts have been made to develop mean field 
models either in r space or in a basis with an improved asymptotic behavior at 
large distances [129, 63, 130, 65, 131, 132, 133]. In particular, an expansion in a 
Woods-Saxon (WS) basis was proposed as a reconciler between the HO basis and the 
integration in coordinate space [65]. Woods-Saxon wave functions have a more realistic 
asymptotic behavior at large r than the harmonic oscillator wave functions. A discrete 
set of Woods-Saxon wave functions is obtained by using box boundary conditions to 
discretize the continuum. It has been shown in Ref. [65] for spherical systems that the 
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solution of the relativistic Hartree equations in a Woods-Saxon basis is equivalent to 
the solution in coordinate space. The Woods-Saxon basis has been used recently for 
both nonrelativistic [134] and relativistic Hartree-Fock-Bogoliubov theories [114] with a 
finite range pairing force. 

Over the past years, lots of efforts have been made to develop a deformed relativistic 
Hartree (RH) theory [135] and a deformed relativistic Hartree-Bogoliubov theory in 
continuum (DRHBc theory) [108]. As a first application, halo phenomena in deformed 
nuclei have been investigated within the DRHBc theory [30]. The detailed theoretical 
framework is available in Ref. [31]. 

For axially deformed nuclei with spatial reflection symmetry, the potentials ¥s(r) 
and Xo(r) in Eq. (14) and various densities can be expanded in terms of the Legendre 
polynomials [136], 


F= X fa(r)Pa(cos0), A= 0,2,4,---, (43) 
with 
i) = “At : I. d(cos 0) f(r) Py(cos 6). (44) 


The quasiparticle wave functions wy, and wy, are expanded in terms of spherical 
Dirac spinors Ynem(rsp) with the eigenvalues cny obtained from the solution of a 
Dirac equation h‘,’ containing spherical potentials S% (r) and V(r) of Woods-Saxon 
shape [65, 128]: 


pra = St Damen PSD) (45) 
Vv (rsp) = b3 Ug a Price PSD): (46) 


NK 
The basis wave function reads 


1 Ganr) (QS) 
ae —Fag(?) Yin (Os) | ea 


where G,,(r)/r and F,,.(r)/r respectively the radial wave functions for the upper 
and lower components. The state Ønsm(rsp) is the time reversal state of Ynam(Tsp). 
The spherical spinor Yngm is characterized by the radial quantum number n, angular 
momentum j, and the parity 7; the latter two are combined to the relativistic quantum 
number x = a(—1)/+!/2(j7 + 1/2) which runs over positive and negative integers 
k = +1,+2,---. For the spinor spherical harmonics )j,, and Vo l = j + $sign(«) 
and 1 = j— ssign(k). These states form a complete spherical and discrete basis in Dirac 
space (see Ref. [65] for details). Because of the axial symmetry the z-component m of 
the angular momentum j is a conserved quantum number and the RHB Hamiltonian 
can be decomposed into blocks characterized by m and parity m. For each (mz)-block, 
solving the RHB equation is equivalent to the diagonalization of the matrix 


( T i ae. A ) Ge = Fe ei (48) 
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where 
ea (aya Vea ye) (49) 
and 
A= (HG eu) = ((n4sme|hp|n'x'mn)), (50) 
B= Cass = ((nama|Alniatmn)) . (51) 


Further details are given in Ref. [31]. 

After solving the RHB or RHFB equations self-consistently, one gets the 
quasiparticle wave functions and meson fields of a nucleus. Then various physical 
quantities can be calculated. Here we take the deformed RHBc as an example and show 
how the total energy, the radius, and the deformation parameters are calculated [31]. 

The total energy of a nucleus is 


E = E nacion F Eg =F Ez T Ep F Ez F Eein, 
2(A — Ex)Ui(m x) Epai 


mT) 
“i J Pr lootjes) +u) 
-5 I a A 
-5 f drg plr) [pe (r) — p3 (r)] 
-5 / dr Agp@(r) + Eem., (52) 


where 
2 
hem | Privo = D (E) (53) 
and U (ø) is the nonlinear self-coupling term of o meson. 
For a zero-range force the pairing field A(r) is local, and the pairing energy is 
calculated as 


1 3 
Bce f dPrr(r)A(r). (54) 


Note that in the DRHBc theory, a smooth cutoff with two parameters, similar to the soft 
cutoff proposed in Ref. [137], has been introduced to regularize the zero-range pairing 
force; for the details, the readers are referred to Refs. [30, 31]. The center of mass 
correction energy 


Eom. = —z7—(P”), (55) 


is calculated after variation with the wave functions of the self-consistent solution [138, 
139] or in the oscillator approximation 


Eom. = -= x 41 x A! MeV, (56) 
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The root mean square (rms) radius is calculated as 


Ress ( J Br oo) 


( f dr [r4p7,_9(r)| ) , (57) 


where T represents the proton, the neutron, or the nucleon. The rms charge radius is 


calculated simply as rg, = r4 + 0.64 fm?. The intrinsic multipole moment is calculated 
by 


= Fal fe [r*p5,a(r)] . (58) 


(59) 


where N, refers to the number of neutron, proton, or nucleon. 


3. Neutron halos and giant halos 


Halo is one of the most interesting exotic nuclear phenomena. There are many new 
features in the halo nuclei because of their large spatial distribution. To give a proper 
theoretical description of the halo phenomena, it is crucial to treat properly the coupling 
between bound states and the continuum due to pairing correlations by either the finite 
range or the zero range pairing forces and the very extended spatial density distributions. 
By assuming spherical symmetry, the RCHB theory [39] provides a proper treatment 
of pairing correlations in the presence of the continuum and an excellent description 
for exotic nuclei [34, 28, 140, 29, 141, 38, 101, 35, 142, 143, 144, 145]. The first 
microscopic self-consistent description of halo in ''Li has been achieved by using the 
RCHB theory [34]. Later, the giant halos in light and medium-heavy nuclei have been 
predicted [28, 39, 29, 144, 145]. The RCHB theory has been generalized to treat the odd 
particle system [101] and extended to including hyperons [104]. Recently, the relativistic 
Hartree-Fock-Bogoliubov (RHFB) theory in continuum has also been developed and 
used to study neutron halos [107, 114]. 

In this Section, progresses on the description of halos in spherical nuclei will be 
reviewed, including the self-consistent description of the halo in 'Li, the giant halo and 
the extension of nuclear landscape, halo and giant halo from the non-relativistic and 
relativistic approaches, and halos from the RHFBc, etc. 


3.1. Neutron halo in “Li 


The ground state properties of Li isotopes have been investigated by using the RCHB 
theory [34]. A satisfactory agreement with experimental values is found for the binding 
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Figure 1. Calculated and experimental density distributions in ''Li and °Li. The 
solid line shows the result of ''Li while the dashed line corresponds to the calculation 
for °Li. The shaded area gives the experimental results with error bars. Taken from 
Ref. [34]. 


energies and the radii of °''Li. The matter radius, showing a considerable increase 
from the nucleus °Li to ''Li, has been reproduced by the RCHB theory. In contrast to 
the earlier mean field calculations, e.g., Refs. [146, 147, 96], the results given in Ref. [34] 
were obtained without any artificial modifications of the potential. 

In Fig. 1 the neutron density distributions of °Li and “Li are shown. It is clearly 
seen that the increase of the matter radius is caused by the large neutron density 
distribution in Li. Its density distribution is fully coincident with the experimental 
density shown with its error bars by the shaded area. 

By examining the mean fields for protons and neutrons and the single particle levels 
in the canonical basis, it has been shown that the halo in 'Li is formed by Cooper- 
pairs scattered in the two levels 1p,/2 and 2581/2 with the former below the Fermi surface 
and the latter in the continuum; the pairing interaction couples those levels below the 
Fermi surface with the continuum. The RCHB theory provides a much more general 
mechanism for halo: a halo appears if the Fermi level is close to the continuum and 
around the Fermi level there are single particle levels with low orbital angular momenta 
and correspondingly low centrifugal barriers. 


3.2. Giant halos and the extension of nuclear landscape 


The halo nuclei observed so far are composed of one or two halo nucleons only. The 
RCHB theory has been used to study the influence of correlations and many-body 
effects in nuclei with a larger number of neutrons distributed in the halo channel. A 
giant neutron halo has been predicted in Zr isotopes close to the neutron drip line [28]. 
It is formed by up to six neutrons scattered to weakly-bound or continuum states beyond 
the !??Zr core with the magic neutron number N = 82. 

Before we proceed, it is worthwhile to discuss more about the characterization 
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of a halo in medium-heavy nuclei which is one of the most relevant topics in the 
study of the halo phenomenon. The critical conditions for the appearance of a halo 
require the existence of two groups of orbitals with significantly different asymptotic 
slopes corresponding to the core and halo nucleons. The one-nucleon or multi-nucleon 
separation energies are close to zero, and there is a significant energy gap for the core 
and valent nucleons in the spectrum. 

Although the theoretical description of light halo nuclei is well under control, as 
discussed in Refs. [148, 149], existing definitions and tools are often too qualitative and 
the associated observables are incomplete for heavier ones. There have been lots of efforts 
in quantifying halos by examining the density profiles or the particles in the classical 
forbidden area obtained from mean field calculations [150, 151, 152]. In Ref. [148], 
a quantitative analysis method is proposed in which the halo part of the density is 
identified as the region beyond a radius rọ where the core density is one order of 
magnitude smaller than the halo one. The method does not require an a priori separation 
of the density into a halo and a core part. Once rọ is fixed, two criteria are introduced to 
characterize halo systems: (i) the average number of fermions participating in the halo 
and (ii) the influence of the latter on the system extension [148]. This method has been 
used to characterize halo features not only in nuclei with the spherical Hartree-Fock- 
Bogoliubov model but also in other finite many-fermion systems like atom-positron/ion- 
positronium complexes [149]. Similar attempts are highly demanded for systematic 
studies of nuclear halos within the frame of both relativistic and non-relativistic DFT. 
In this Review, we mainly analyze the halo properties based on a combined analysis of 
the following features: i) the deviation of the rms radius from the roN1/? law; ii) the 
large spatial extension of the density profile; and iii) the occupation of weakly bound 
and/or continuum orbitals in the canonical basis. For those weakly bound or continuum 
orbitals, we also discuss in detail their contribution to the halo. 

In the upper panel of Fig. 2, the rms radii of the protons and neutrons are presented 
for the Zr isotopes. There is a clear kink for the neutron rms radius at the magic neutron 
number N = 82. This kink can be understood by examining the single particle levels 
in the canonical basis which is shown in the lower panel of Fig. 2. Going from N = 70 
to N = 100, there is a big gap above the 1h11/2 orbital. For N > 82, the neutrons 
are filled to the levels in the continuum or weakly bound states in the order of 3p3/2, 
2 f7/2, 3P1/2, 2fs5/2 and 1hg/2. The neutron chemical potential is given by a dashed line. 
It approaches rapidly the continuum threshold already shortly after the magic neutron 
number N = 82 and crosses the continuum at N = 100 for the nucleus to “°Zr. In 
this region the chemical potential keeps a very small but negative value for a large A 
range. This means that the additional neutrons are added with a very small, nearly 
vanishing binding energy at the edge of the continuum provided by pairing correlations. 
The total binding energies E for the isotopes above !??Zr are therefore almost identical. 
This has been recognized in Ref. [95]. However, the RMF calculations with the BCS 
approximation expanded on an oscillator basis is definitely not appropriate for exotic 
nuclei close to the continuum limit. 
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Figure 2. Upper part: Root mean square radii for neutrons and protons in Zr isotopes 
close to the neutron drip line as a function of the mass number A. Lower part: single 
particle energies for neutrons in the canonical basis. The dashed line indicates the 
chemical potential. Taken from Ref. [28]. 


From the occupation probabilities in the canonical basis, one can get the number 
of halo neutrons. It was found that there are 2 valence neutrons in '**Zr, 4 in !°Zr, 
6 in Zr, ..., roughly 5 in 8Zr, and roughly 6 in “°Zr where the neutron drip line 
is reached. Note that this does not mean that there are around six neutrons in the 
halo region of the nuclear density in °Zr. We emphasize that these figures, 2, 4, 5, or 
6, are numbers of “valence” neutrons, just as it is well recognized that there are two 
“valence” neutrons in “Li. With the very large neutron rms radii of these systems, one 
can estimate the number of valence neutrons which fill in the same volume outside the 
12277 core if packed with normal neutron density. This number is 24 for !°*Zr and 34 
for 4°Zr. This phenomenon is therefore clearly a neutron halo and was called a Giant 
Halo because of the large number of particles participating in the halo formation [28]. 

It is very difficult for experimentalists to study giant halos in Zr isotopes because 
they are too heavy to be synthesized by the RIB facilities at present. Thus it is useful 
to investigate the giant halo phenomena in lighter nuclei which are easily accessible 
with available facilities. In Refs. [29, 145], ground state properties of all the even-even 
O, Ca, Ni, Zr, Sn, and Pb isotopes ranging from the proton drip line to the neutron 
drip line were investigated by using the RCHB theory with the effective interaction 
NLSH [153]. Recently, a systematic calculation for nuclei above oxygen has also been 
performed with the effective interaction PC-PK1 [118], aiming to investigate the global 
impact of the continuum for the nuclear boundary; the results ranging from O to Ti 
isotopes have been published in Ref. [5]. In such a systematic study, it is found that 
the appearance of giant halos, the contribution of continuum and boundaries of nuclear 
chart are tightly connected. We note that the functional PC-PK1 has proved to be very 
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successful in describing the isospin dependence of the nuclear masses [154], the Coulomb 
displacement energies between mirror nuclei [155], fission barriers [156, 157] and nuclear 
rotations [92, 93, 91, 112], ete. 

For nuclei far from the valley of stability and with small nucleon separation energy, 
the Fermi surface is very close to the continuum threshold and the valence nucleons 
may extend over quite a wide space to form low density nuclear matter. Therefore 
the nucleon separation energies are sensitive quantities to test theory and examine the 
halo phenomena. The two-neutron separation energies for O, Ca, Ni, Zr, Sn, and Pb 
isotopes have been investigated to examine their connection with systematic behavior 
of the nuclear size [29, 144, 145]. 


r, [fm] 


0 20 40 60 80 100 120 140 160 
N 


Figure 3. (Color online) (a) The two-neutron separation energies Sən calculated by 
RCHB theory [39] with NLSH [153] (Open symbols) in comparison with data available 
(solid ones) for the even O, Ca, Ni, Zr, Sn, and Pb isotopes against the neutron number 
N. (b) The root mean square (rms) neutron radii r» for even O, Ca, Ni, Zr, Sn, and 
Pb isotopes by RCHB calculations against the neutron number N. The curve r9N1/3 
with ro = 1.139 fm has been included to guide the eye. Taken from Ref. [29] with 
modified color. 


In Fig. 3, similar to Ref. [29], the two-neutron separation energies Sən and the 
neutron radii r„ from the RCHB calculation for the even-even nuclei of the O, Ca, Ni, 
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Zr, Sn, and Pb isotope chains have been shown. The predicted r, curve using the 
simple empirical equation r,, = ro N'/3 with ro = 1.139 fm normalizing to 2°°Pb is also 
represented in the figure. This simple formula for r, agrees with the calculated neutron 
radii with exceptions in neutron rich Ca and Zr isotopes. 

The comparison between experimental and calculated So, is shown in Fig. 3. The 
experimental magic or submagic numbers N = 20,28, and 40 are reproduced. The Son 
values for several nuclei in exotic Ca and Zr isotopes are extremely close to zero. As 
discussed in Refs. [28, 29], if taking Ca or '”*Zr respectively as a core, the valence 
neutrons will gradually occupy the loosely bound states and the continuum above the 
sub-shell of N = 40 in Ca or N = 82 in Zr isotopes. Such typical behavior of Sən can 
be taken as an evidence of the occurrence of giant halos in Ca chain and Zr chain [28]. 

In Fig. 3, it is very interesting to see that r, follows the M! systematics well 
for stable nuclei although their proton numbers are quite different. Near the drip line, 
distinct abnormal behaviors appear at N = 40 in Ca isotopes and at N = 82 in Zr 
isotopes. As pointed in Ref. [29], nuclei having the abnormal r,, increase correspond 
to those having small Ssn, which provides another evidence to the emergence of giant 
halos. The increase of r, in exotic Ni, Sn and Pb nuclei is not as rapid as those in 
Ca and Zr chains. The regions of the abnormal increases of the neutron radii are just 
the same as those for Sən. Both behaviors are connected with the formation of giant 
halo [28, 29, 144, 145]. 


3.8. Halos from relativistic and non-relativistic approaches 


By using the RCHB theory, predictions of halo in the Ne [88, 158], Na [101, 158], Ca 
[29] and Zr [28] isotopes near the neutron drip line have been made. In particular, in 
Ca and Zr isotopes near the neutron drip line, the halo phenomena were addressed as 
giant halo [28, 29]. In the non-relativistic HFB approaches, Skyrme calculations with 
the parameter set SLy4 predicted the halo in Sn and Ni [152]. It was also claimed that 
the pairing gaps have an effect to reduce the halo at the neutron drip line [152, 159, 160). 

In Refs. [161, 162], the phenomena of the giant halo and halo of the neutron-rich 
even-Ca isotopes have been investigated and compared between the frameworks of the 
RCHB and the Skyrme HFB calculations. With the two parameter sets for each of the 
RCHB and the Skyrme HFB calculations, it has been found that although the halo 
phenomena exist for Ca isotopes near the neutron drip line in both calculations, the 
halo from the Skyrme HFB calculations (SkM*) starts at a more neutron-rich nucleus 
than that from the RCHB calculations, and SLy4 does not predict a halo, because the 
drip line is much closer to the stable region than the other parameter sets. 

Figure 4 shows systematics of the rms radii of neutrons 7,, and protons 7, of even- 
Ca from N = 30 to the two-neutron drip lines. The increase in the curvature of 7, 
indicates the halo. For Ca, the number of nucleons in continuum, Nj;,, of ©?~7?Ca is 0.6- 
2.2, of which the average is 1.7, in the NL-SH calculation [29]; the corresponding value 
of the SkM* calculation is always smaller than 0.5 (the average 0.27). As there are more 
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Figure 4. The neutron and proton rms radii 7, and 7, of even-Ca isotopes in RCHB 
calculations with NLSH and PK1 and the Skyrme HFB calculations with SkM* and 
SLy4 in the neutron-rich region. Taken from Ref. [161]. 


than two halo neutrons, it was referred to as giant halo [28, 29]. The SkM* calculation 
predicts the halo from N = 52, and the RCHB calculations show the gradual occurrence 
of the giant halo. Apparently, the starting nucleus of the halo of SkM* corresponds to 
the extra lowering of the Sən. The SLy4 calculation does not predict the halo, because 
the particle-stable region ends at N = 46. The large difference between 7, and 7, shown 
in Fig. 4 was identified as neutron skin in the region with no halo. 

The RCHB calculations have larger neutron rms radii systematically in N > 40 
than the Skyrme HFB calculations. This difference comes partially from the occupation 
of the neutron 3581/2 orbital, which causes 50% of the maximum difference in the neutron 
rms radii among the four calculations at ®°Ca. The difference in the neutron 351/2 orbital 
is shown in Fig. 5, which illustrates the radial wave functions squared. The tail of 351/2 
from PK1 has appreciably extended distribution than the others, and as a consequence, 
the amplitude of PK1 is smaller in the inner region. The number and locations of the 
nodes are the same for the four curves. 

As shown in the previous Subsection and in Ref. [28], the occupation of 3p;/2 and 
3p3/2 orbitals are responsible for the formation of giant halos in Zr isotopes. In Ref. [163], 
the Hartree-Fock-Bogoliubov (HFB) approach with Skyrme interactions SLy4 and SkI4 
has been used to study halos in Ca and Zr isotopes. It was found that the appearance 
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neutron 3s1/2 


Figure 5. Radial wave function squared |7,|? (canonical basis) of the neutron 3s1/2 
of ©°Ca. They are normalized as f dr |w,|? = 1. Taken from Ref. [161]. 


of giant halos depends sensitively on the effective interaction adopted: SkI4 predicts a 
neutron halo in the Zr chain with A > 122 due to the weakly bound orbitals 3p;/2 and 
3P3/2- 

In the above mentioned work, the box boundary condition is adopted and the 
continuum is discretized. Although the convergence of the halo properties with respect 
to the box size has always been achieved, it is still desirable to examine the influence 
of the asymptotic behavior of the radial wave functions of those valence orbitals on the 
formation of halo or giant halo. The Skyrme HFB equations in r space have been solved 
for spherical systems with correct asymptotic boundary conditions for the continuous 
spectrum [164] or the Green’s function technique [55, 57]. In Ref. [164], it was shown that 
close to the drip line the amount of pairing correlations depends on how the continuum 
coupling is treated. In Ref. [55], it was found that, in even-even N = 86 isotones in 
the Mo-Sn region, the l = 1 broad quasiparticle resonances persist in feeling the pairing 
potential and contribute to the pairing correlation even when the widths are comparable 
with the resonance energy. In Ref. [56], the investigation of the pair correlation in the 
tail of the giant halo revealed that the asymptotic exponential behavior of the neutron 
pair condensate is dominated by nonresonant continuum states with low asymptotic 
kinetic energy. 


3.4. Halos with relativistic Hartree-Fock-Bogoliubov theory in continuum 


In spite of successes of RCHB theory, there are still a number of questions needed to 
be answered in the conventional CDFT, e.g., the contributions due to the exchange 
(Fock) terms and the pseudo-vector 7-meson. Earlier attempts to include the exchange 
terms with relativistic Hartree-Fock (RHF) method led to underbound nuclei due to the 
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missing of the meson self-interactions [165, 126]. For a long time, the RHF theory failed 
in a quantitative description of nuclear systems [166, 167, 168, 169, 170, 126, 171, 172). 
With the development of the computational facilities and numerical techniques, the 
density dependent relativistic Hartree-Fock (DDRHF) theory has been developed and 
shown significant improvements in a quantitative description of nuclear phenomena 
[127, 78, 113, 173, 174, 175, 176, 177] with comparable accuracy. With a number of 
adjustable parameters comparable to those of RMF Lagrangians, the DDRHF theory 
can give a equally good description of nuclear systems without dropping the Fock 
terms. Furthermore, important features like the behavior of neutron and proton effective 
masses |178] can be interpreted. 

In DDRHF theory [127], the tensor forces due to 7 and p meson exchanges can 
be naturally taken into account and have brought significant improvements on the 
consistent description of shell evolution [173, 179, 180, 181]. The relativistic Hartree- 
Fock-Bogoliubov (RHFB) theory in continuum with density dependent meson-nucleon 
couplings is developed in Refs. [107, 114]. In the descriptions of nuclear halos, by using 
the finite range Gogny force D1S in the pairing channel, systematic RHF Be calculations 
are performed both for stable and weakly bound nuclei. It is demonstrated that an 
appropriate description of both mean field and pairing effects can be obtained within 
RHFBc theory [107, 114]. 

In Fig. 6, the single-neutron separation energies Sn of Sn isotopes from ''Sn to 
'88Sn (left panels) and the single-proton separation energies S, of N = 82 isotones from 
130Cd to 13Lu (right panels) are given. The odd-even differences on the single-nucleon 
separation energies reflect the effects of the pairing correlations. In Fig. 6, RHFBc with 
PKA1 [113] and PKO1 [127] as well as RHB with DD-ME2 [182] present comparable and 
satisfied quantitative agreements with the data for both isotopes and isotones. From 
Fig. 6, one can find some systematics in the results. On the neutron rich side, i.e., 
after '°’Sn for Sn isotopes and before '3¢Gd for N = 82 isotones, PKA1 shows better 
agreements than PKO1 and DD-ME2. This is mainly due to the improvements in the 
shell structures with exchange terms included, especially the shell closures beyond 50 
and 82, i.e., 58 and 92, which may change the level density and then influence the pairing 
effects. On the proton rich side, these three approaches provide comparable agreement 
with the data. 

The Ce isotopes have been chosen to demonstrate the nuclear halo phenomenon 
within the RHFBc theory [107, 114]. In Fig. 7 the nuclear matter distributions (left 
panels) and neutron canonical single particle configurations (right panels) for Ce isotopes 
close to the drip line are shown. The neutron drip line calculated here with PKA1 [113] 
is N = 140 whereas the calculations with PKO1 [127] and DD-ME2 [182] predict a 
shorter one as N = 126. As shown in Fig. 7(a), the neutron densities become more and 
more diffuse after the isotope Ce (N = 128), a direct and distinct evidence of halo 
occurrence. 

From Fig. 7(b) one can see that such extremely extensive matter distribution, e.g. 
in '*8Ce, is mainly due to the low-/ states, namely the halo orbitals 1451/2, v3ds/2 and 
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Figure 6. (Color online) Single-nucleon separation energies along Sn isotopic (Sn: left 
panels) and N = 82 isotonic (Sp: right panels) chains. The results are calculated by 
RHFB with PKA1 [113], PKO1 [127], and by RHB with DD-ME2 [182], in comparison 
to the experimental data [183]. Taken from Ref. [114]. 


v3d3/2. From the occupations of the halo orbitals Nhao in Fig. 7(d), the halos were 
predicted in °Ce, 188Ce and !°°Ce and giant halos may appear in !°?Ce, 194Ce, 1%Ce 
and !°8Ce because more than two neutrons are occupying the halo orbitals. Similar 
conclusions can be obtained from the neutron numbers lying beyond the sphere with 
the radius r = 10 fm [N;s10 fm in Fig. 7(d)], which is large enough (the neutron matter 
radius r„ = 6.2 fm in '*8Ce) for halos. Even extending to r = 16 fm, which is sometime 
taken as the radial cut-off in the calculations of stable nuclei, there are still some amount 
of neutrons lying beyond this sphere in the isotopes from '?Ce to !°Ce. 

In Fig. 7(c), one finds that the halo orbitals (v481/2, v3d5/2 and v3d3/2) are located 
around the particle continuum threshold where they are gradually occupied. For the 
isotopes beyond '*4Ce (N = 126), the Fermi levels (in open circles) approach the 
continuum threshold rather closely such that the stability of these halo isotopes becomes 
sensitive to pairing effects. Nearby the low-/ states, there are the high-I states v2g9/2 
and v297/2. Because of the relatively large centrifugal barrier for g-orbitals they do not 
contribute much to the diffuse neutron distributions. Nevertheless, the existence of the 
high-/ states nearby halo orbitals is still particulary significant, because it leads to a 
rather high level density around the Fermi surface, and evidently the pairing effects are 
enhanced to stabilize the halo isotopes. 

Evidence for the existence of a halo can also be investigated by the systematic 
behavior of nuclear bulk properties such as radii. In Fig. 8, the isospin dependence of 
the neutron skin thickness (r, — rp) calculated in density dependent RHFB theory in 
continuum using PKA1 for Ca, Ni, Zr, Sn, and Ce isotopes are given [107]. Taking the 
stable nuclei as references (shown as the dashed lines), continuously growing deviations 
are found in the Ca, Zr, and Ce chains until the neutron drip line. This can be considered 
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Figure 7. (Color online) (a) Neutron and proton densities, (b) the relative 
contributions of different orbitals to the full neutron density in !°8Ce, (c) neutron 
canonical single particle energies occupation probability (in x-error bars) and Fermi 
energy Ep (in open circles), and (d) neutron numbers filling in the halo orbitals 
451/2, 3d5/2, and 3d3/2 (Nhalo), and the ones lying beyond the spheres with the 
radii r = 10, 11, 12, 13, 14, 15, 16 fm, respectively N;>10 fm, Nr>11 fm, Nr>12 tm; 
Nr>13 fm, Nr>14 fm, Nrsis fm and Nr>16 tm. The results are calculated by RHFB 
with PKA1 [113] plus the Gogny pairing force D1S [124]. The spherical box radius is 
adopted as R = 28 fm. Taken from Ref. [107]. 


as evidence of a halo. Despite the deviations in the mid-region, both Ni and Sn nuclei 
show a similar isospin dependence in both stable and exotic nuclei regions, which indicate 
a neutron skin since the growth of a halo is interrupted. Compared with Ni, the Sn 
isotopes show even a much weaker skin effect. 

The systematics for the Ca and Zr chains from the RHFBc calculations are similar 
to those from the RCHB calculations [29]. For the Ce isotopes, RHFBc calculations 
with PKA1 show clear evidence for the existence of halo structures. The calculations of 
RHFBc with PKO1 [127] and RCHB with DD-ME2 [182] predict that the isotopic chain 
ends at N = 126, before the halo occurrence as predicted by RHFB with PKA1. This 
deviation is interpreted by the shell structure evolution in Fig. 7(c), where much weaker 
shell are predicted with PKA1 in comparison with PKO1 and DD-ME2. As shown in 
Fig. 7(c), the neutron shell gap (N = 126) between 1113/2 and v2g9/2 states is close to 
the particle continuum threshold, which might essentially influence the stability of the 
drip line isotopes. 

From Ref. [107], one may conclude that the Fock term itself has little direct impact 
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Figure 8. (Color online) Neutron skin thickness r„ — rp for Ca, Zr, Ni, Sn, and Ce 
as a function of the isospin (N — Z). The r, and rp are respectively the neutron and 
proton root mean square radii. The results are calculated by RHFB with PKA1 [113] 
plus the paring force D1S [124]. Taken from Ref. [107]. 


on the formation of a nuclear halo. However, the Fock term may influence the shell 
structure, particularly those orbitals close to the Fermi surface in drip-line nuclei, thus 
affecting indirectly the formation of a potential halo. 


3.5. Halos in exotic hyper nuclei 


Motivated by knowledge of A-N interaction and understanding on giant halo [28, 29, 145], 
it is interesting to study the possible appearances of halos in exotic hyper nuclei. 

Two neutron separation energies Sən for normal nuclei, single-A hyper nuclei and 
double-A hyper nuclei of Ca isotopes, labelled by A = 0, A = 1 and A = 2, respectively, 
from the proton drip line to neutron drip line were investigated [104]. It was found that 
one or two A hyperons lower the Fermi level due to the attractive A-N potential but 
keep the neutron shell structure unchanged. Therefore, the neutron drip line is pushed 
outside from N = 52 in Ca isotope chain to N = 54 in hyper isotope chain. Meanwhile, 
giant halos due to pairing correlation and the contribution from the continuum still exist 
in Ca hyper nuclei similar to that in Ca isotopes [29]. This is a slight but rewarding 
step for exploring the limit of drip line nuclei on the basis of the giant halo. For more 
details, see Ref. [104]. 

Apart from neutron halos in hyper nuclei, as A hyperon is less bound than the 
corresponding nucleon in nuclei, it is worth investigating the existence of hyperon halos. 

In Ref. [184, 185], a A halo was tentatively suggested in a RMF plus BCS model. 
However, since the normal BCS method will have unphysical solution by involving 
baryon gas [39], such prediction needs further study. In Ref. [103], hyper carbon isotopes 
are studied by the RCHB theory. By investigating the baryon density distributions in 
hyper nuclei °C, 14C, and 3iC, it was found that with up to two A hyperons added 
to the core °C, the nucleon density distributions remain the same and hyperon density 
distributions at the tail are comparable with those of the nucleons. An intriguing 
phenomenon appears in 3,C where the hyperon density distribution has a long tail 
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extended far outside of its core 14C. This is a signature of hyperon halo due to the 
weakly bound state 1p3/9 in 44C which has a small hyperon separation energy and a 
density distribution with long tail [103]. 


4. Surface diffuseness 


From the mean field point of view, the properties of nucleons in an atomic nucleus 
are determined by the mean potential which stems from their interaction with the other 
nucleons and is determined by nucleon densities in the self-consistent models. Therefore, 
the study of the isospin dependence of the mean potential and the density distribution, 
both become highly diffuse near the particle drip line, is crucial to understand exotic 
nuclei. For finite nuclei, the diffuseness of the nuclear surface provides a measure of the 
thickness of the surface region. The surface diffuseness is intimately related to the spin- 
orbit splitting [142], the evolution of single particle shell structure [28], the pseudospin 
symmetry [76, 77, 186] and spin symmetry [83, 186], and the nuclear surface energy and 
thus nuclear masses [187]. In this Section, we will review recent progresses on nuclear 
surface diffuseness and the corresponding impact on nuclear shell structure and size. 


4.1. Spin-orbit splitting and shell structure 


The RCHB theory has been used to systematically study the isospin dependence of the 
mean potential for the Sn isotopes from proton drip line to neutron drip line and the 
spin-orbit splitting for the whole isotope chain was examined in details [142]. As an 
example, the spin-orbit splitting 


E Eyjat—-1/2 — Pij=1+1/2 


Pis 21+ 1 (60) 
versus the binding energy: 
1+ 1)Fyai-12 + lE- 
E= ( )Eij=1-1/2 lj=l+1/2 61) 


2i+1 
in “Sn with A = 110,120,---,170 are given in Fig. 9 for the neutron spin-orbit 
partners (1d3/2,1d5/2), (1g7/2,1g9/2), (Ati1y2, 1413/2), (1p1/2, 1p3/2), (fsy2, 1 f7/2) and 
(1hg/2, Lh11/2), and the proton spin-orbit partners (1d3/2, 1d5/2) and (1f5/2, 1 f7/2). 

It is very interesting to see that the spin-orbit splittings for the neutron and proton 
are very close to each other, at least for (1d3/2,1d5/2) and (1fs/2, 1fz/2) cases. The 
splitting decreases monotonically from the proton drip line to the neutron drip line. 
To understand its mechanism, it is very helpful to examine the origin of the spin-orbit 
splitting in the Dirac equation. 

For the Dirac nucleon moving in scalar and vector potentials, its equation of motion 
could be de-coupled and reduced for the upper component and the lower component, 
respectively. For the spin-orbit splitting, the Dirac equation can be reduced for the 
upper component as following: 

d? 1 d(2M — Xo + £s) d 


Poe o UT 2 
Ge ETM SoBe dr ap G8) 
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Figure 9. The neutron spin-orbit splitting Ers = (Eyj=i—1/2 — Exjai4i/2)/(2! + 1) 
versus the mass number A in Sn isotopes for neutron (1d3/2,1ds5/2), (197/2, 199/2), 
(li11/2, 1¢13/2), (p12, 1p3/2), (1fs/2,1f7/2) and (1hg/2,1h11/2) orbitals and proton 
(1d3/2, 1ds/2) and (1fs/2,1f7/2) orbitals, respectively. Taken from Ref. [142]. 


7 pags, 7 1 PULM = Ag 
p2 E+2M—%)+NUgr dr í 
= — (E +2M — Xo + X5)(E — Xo — X5)GË (r). (62) 
The spin-orbit splitting is due to the corresponding spin-orbit potential 
1 k d(2M — Xp + Xs) (63) 
E+2M — Xo+ Xsr dr 


with some proper normalization factor. The spin-orbit splitting is energy dependent 
and depends on the derivative of the potential 2M — Xo + Ug. Therefore the so-called 


c d(2M — Yo +È 
spin-orbit potential, Vj, = KAUM T 
Fi 


Tr , is introduced. 

The derivatives of the neutron potentials No(r) — Ug(r) for “Sn, “°Sn, and tSn, 
are given in the upper panel of Fig. 10. For both proton and neutron, the derivatives 
of the potentials Uo(r) — Us(r) are almost the same. As the potential Uo — Us is a 
big quantity (~ 700 MeV), the isospin dependence in the spin-orbit potential could be 
neglected. That is the reason why the spin-orbit splitting for the neutron and proton is 
very close to each other in Fig. 9. From ''°Sn to tSn, the amplitude of the derivative 
for Xo — Ug decreases monotonically due to the surface diffuseness. 

The decline of the spin-orbit splitting is from the diffuseness of the potential or 
the outwards tendency of the potential. The neutron and proton potentials Uo + Ug 
are given in Fig. 11 for “Sn with A = 110,120,---,170. It is seen that the depth of 
the neutron potential decreases monotonically from the proton drip line to the neutron 
drip line and the surface moves outwards. The inserted figures give the radii Ro at 
which Mop + Ns = —10 MeV as functions of the mass number. With the increase of 
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Figure 10. Upper: The derivative of the neutron potentials Xo(r) — Ss(r) for 4Sn 
with A = 110, 120,---,170; Lower: The neutron potentials Xo(r)— ¥s(r) for 4Sn with 
A = 110,120,---,170. In order to examine the surface diffuseness more clearly, the 
radii Ro at which “o(Ro) — Us(Ro) = 100 MeV has been given as an inserted figure. 
Taken from Ref. [142]. 


the neutron number, the depth of the proton potential increases monotonically but is 
pushed towards outside as well due to the proton-neutron interaction. 

Similar neutron potentials Ng — Ug for 4Sn with A = 110, 120,---,170 are given in 
the lower panel of Fig. 10 and its inserted figure gives the radii Rọ at which Xọ—¥s = 100 
MeV as a function of the mass number. As seen above, the spin-orbit splitting is related 
with the derivative of the potential No — “is. The surface diffuseness happens for both 
the vector and scalar potential as shown in both My — Ng and Np + Ng. 

Similar RCHB calculations have been carried out for all the nuclei in Na isotopes 
with mass number A ranging from 17 to 45 with effective interaction NLSH [28]. Apart 
from the diffuseness of the potential, the shell evolution with neutron number has been 
examined there as well. In Fig. 12, the microscopic structure of the single particle 
energies in the canonical basis [188, 39] is given. In the left panel of Fig. 12, the single 
particle levels in the canonical basis for the Na isotopes with an even neutron number are 
shown. The level density becomes dense with the neutron number due to the decrease 
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Figure 11. Left: The neutron potentials Xo(r) + Es(r) for 4Sn with A = 
110, 120,---,170. In order to examine the surface diffuseness more clearly, the radii 
Ro at which Xo(Ro) + Es(Ro) = —10 MeV has been given as an inserted figure; 


Right: The proton potentials Xo(r) + Ss(r) for 4Sn with A = 110,120,---,170. 
In order to examine the surface diffuseness more clearly, the radii Rọ at which 
Xo(Ro) + Us(Ro) = —10 MeV has been given as an inserted figure. Taken from 
Ref. [142]. 


of the spin-orbit potential or diffuseness of the potential. Going from A = 19 to A = 45, 
it can be observed a big gap above the N = 8, N = 20 major shell, and N = 14 sub- 
shell. The N = 28 shell for stable nuclei fails to appear, as the 2p3/2 and 2p;/2 come so 
close to 1f7/2. When N > 20, the neutrons are filled to the levels in the continuum or 
weakly bound states in the order of 1 f7/2, 2p3/2, 2p1/2, and 1fs/2. In the right part, the 
occupation probabilities in the canonical basis of all the neutron levels below Æ = 10 
MeV have been given for Na to show how the levels are filled in nuclei near the drip 
line. The importance of careful treatment of the pairing correlation and the scattering 
of particle pairs to higher lying levels are noted in the figure. 


4.2. Interaction cross section 


The neutron and proton distributions and radii are fundamental quantities to define the 
halo phenomena. However, the radius for unstable nucleus is normally extracted from 
the interaction cross sections by assuming an empirical form (e.g., a Gaussian function) 
for the density distribution [24]. It would be better to use the density distribution 
from the self-consistent models and calculate the interaction cross sections directly from 
Glauber model. 

Such systematic calculations have been carried out for all the nuclei in Na isotopes 
with mass number A ranging from 17 to 45 with the RCHB theory and the effective 
interaction NLSH [101]. The calculated binding energies Æ and the interaction cross 
sections with the Glauber model are shown in Fig. 13. The calculated binding energies 
E are compared with the empirical values [189]. The proton (neutron) drip line nucleus 
has been predicted to be !°Na (Na). The difference between the calculations and 
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Figure 12. Left part: Single particle energies for neutrons in the canonical basis as a 
function of the mass number for Na isotopes. The dashed line indicates the chemical 
potential. Right part: The occupation probabilities in the canonical basis for *°Na. 
Taken from Ref. [101]. 


the empirical values for the stable isotopes is from the deformation, which has been 
neglected in Ref. [101]. 

As the interaction cross sections gz is well reproduced, the density distributions of 
the whole isotopes from the RCHB theory have been examined and the relation between 
the development of halo and shell effect has been studied in detail [101]. It was found 
that the neutron N = 28 shell closure fades due to the lowering of the 2p orbitals. 
The tail in the nuclear density profile depends strongly on the shell structure and the 
development of a halo was connected with the changes in the occupation in the next 
shell or sub-shell close to the continuum limit [101]. 

Similar calculations have been performed for C, N, O and F isotopes up to the 
neutron drip line by the self-consistent RCHB theory [143]. The agreement between the 
calculated one neutron separation energies S» and the available experimental values are 
similar as Na isotopes. A Glauber model calculation for the total charge-changing cross 
section has been carried out with the density distribution obtained from the RCHB 
theory. The measured cross sections for studied nuclei with 1?C as a target are nicely 
reproduced. An important conclusion was found that, contrary to the usual impression, 
the proton density distribution is less sensitive to the neutron number along the isotope 
chain. Instead it is almost unchanged from stability to the neutron drip line. Such 
calculations are quite useful in extracting both the proton and neutron distributions 
inside the nucleus and in defining the neutron/proton skin or halo. 
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Figure 13. Upper part: The interaction cross sections øz of “Na isotopes on a 
carbon target at 950A MeV: the open circles are the result of RCHB calculation and 
the available experimental data (A =20—23, 25-32) are given by solid circles with their 
error bars. The dashed line is a simple extrapolation based on the RCHB calculation 
for 2°~3'Na. Lower part: Binding energies for Na isotopes, the convention is the same 
as the upper part, but the RCHB result for particle unstable isotopes are indicated by 
triangle. Taken from Ref. [101]. 


4.3. Surface diffuseness correction in global mass formula 


For finite nuclei, the diffuseness of nuclear surface is an important degree-of-freedom 
in the calculations of nuclear masses [187]. The density distributions of most stable 
nuclei are of the “neutron skin-type”, with a typical value around 0.5 fm for the surface 
diffuseness. For nuclei near the neutron drip line, halos and giant halos may develop and 
the neutrons may distribute more extensive spatially than the protons, which implies the 
enhanced neutron surface diffuseness for these extremely neutron-rich nuclei. However, 
in global mass calculations, the surface diffuseness of exotic nuclei near the drip lines 
has not been properly considered. Recently, the surface diffuseness from the RCHB 
theory and its impact on the symmetry energy and shell correction have been studied in 
Ref. [187] and the accuracy of the global mass formula has been considerably improved. 

Inspired by the Skyrme energy-density functional, a macroscopic-microscopic mass 
formula, the Weizsacker-Skyrme (WS) formula [190, 191, 192], was proposed with a rms 
deviation around 336 keV with respect to the 2149 measured masses [183] in 2003 Atomic 
Mass Evaluation (AME). The WS models provide the best accuracy for nuclear masses 
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Figure 14. (Color online) (a) Difference between measured and the WS4 calculated 
masses. (b) Difference between WS3 and WS4 calculated masses. The squares and 
crosses denote the nuclei first appearing in AME2012 and the nuclei with |I—Jo| > 0.1, 
respectively. The smooth and the zigzag curves denote the -stability line from the 
Green’s formula and the neutron drip line from the WS4 formula, respectively. Taken 
from Ref. [187]. 


in several mass regions [193, 194] and has been used in r-process simulations [195, 196]. 
In the WS formula, an axially deformed Woods-Saxon potential, with a contant 
surface diffuseness parameter a for all nuclei, was used to obtain the single particle 
levels of nuclei. In Ref. [187], by taking into account the surface diffuseness correction 
for unstable nuclei, a new global mass formula, WS4, has been proposed. The surface 
diffuseness a of the Woods-Saxon potential was given by a = ao (1+ 2e6,), with 
e = (I — Ip)? — I* the correction factor for surface diffuseness, ao the diffuseness of the 
Woods-Saxon potential, Tọ = 0.44/(A + 200) the isospin asymmetry of the nuclei along 
the -stability line described by Green’s formula, and 6, = 1 for neutrons (protons) in 
the nuclei with 7 > Ig (I < Io) and 6, = 0 for other cases. This means that the surface 
diffuseness of neutron distribution is larger than that of protons at the neutron-rich 
side and smaller than that of protons at the proton-rich side. The rms deviation with 
the available mass data [197] falls to 298 keV, crossing for the first time the 0.3 MeV 
accuracy threshold for mass formulas (models) within the mean field framework. 

In Fig. 14(a), it was shown the deviations of the calculated masses from the 
experimental values. For all the 2353 nuclei with Z and N > 8, the deviation is within 
1.23 MeV. In Fig. 14(b), the difference between the WS3 (constant surface diffuseness ag) 
and WS4 calculated masses. For most nuclei, the results of WS3 and WS4 are consistent 
in general (with deviations smaller than one MeV). For nuclei near the neutron drip line, 
the masses given by WS4 are larger than the results of WS3 by several MeV. This is due 
to the enhancement of the nuclear symmetry energy coming from the surface diffuseness 
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effect in the extremely neutron-rich nuclei. 

It is worthwhile to note that with an accuracy of 258 keV for all the available 
neutron separation energies and of 237 keV for the alpha-decay @-values of superheavy 
nuclei, the proposed mass formula will be important not only for the reliable description 
of the r-process of nucleosynthesis but also for the study of the synthesis of superheavy 
nuclei. 


5. Pairing correlations and nuclear size 


In Ref. [159], it is shown that the neutron ground state Hartree-Fock (HF) asymptotic 
density equals 


p(r) x exp (—2ur)/r?, (64) 
with u = /—2me;/h and €x the single particle energy of the least bound / = 0 neutron. 
The corresponding mean square radius deduced from the asymptotic solution is 
h2 


——— 65 
A (65) 


(r?) HF 


which diverges in the limit €% — 0. 
In the presence of pairing, the mean square radius deduced from the asymptotic 

Hartree-Fock-Bogoliubov (HFB) density is 
h? 


with the lowest discrete quasiparticle energy Ep = \/ (ex — A)? + AZ. If the pairing 
gap A, is finite, the radius will never diverge in the limit of small separation energy 
ERS ÀA >O. 

The asymptotic HF and HFB densities characterized by l = 0 orbitals were 
compared in Ref. [159] and it has been emphasized that pairing correlations reduce 


(r?) HEB & (66) 


the nuclear size and an extreme halo with infinite radius cannot be formed in superfluid 
nuclear systems. Here pairing correlations act against the formation of an infinite radius. 
Therefore this effect was called “pairing anti-halo effect” in Ref. [159]. 

In fact, the mechanism leading to a halo in Eq. (65), i.e., the valence nucleons 
occupy the weakly-bound orbital below and very close to the continuum threshold, has 
been used in early interpretation for nuclear halo phenomena [146, 147, 96]. The limiting 
condition €% — 0 and the radius deduced from this l = 0 orbital alone correspond to an 
extremely ideal situation, which is difficult to be found in real nuclei. 

In Ref. [198], the influences of pairing correlations on the nuclear size and on the 
development of a nuclear halo were studied in details by the self-consistent RCHB 
theory [39]. In order to simplify the problem, neutron-rich nuclei with the neutron 
Fermi surface below, between and above two weakly-bound 2p levels are investigated 
with a fixed Wood-Saxon potential, which is followed by the self-consistent calculation 
on the well studied neutron halos in Li and **Ne with low-/ orbitals in the continuum. 
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5.1. Fermi surface and nuclear size 


In order to investigate the impact of nuclear binding on nuclear size for varying 
pairing potential, a spherical Woods-Saxon shape mean field potential fitted to the 
corresponding neutron potentials resulting from a self-consistent calculation for Mg is 
adopted [198]. 
The total neutron radius of the nucleus is given by 
Ry = L S72 + 1) nes ay (67) 
nlj 

which is determined by the rms radii of the orbitals (r),,; and the corresponding 


occupation probabilities U2: 
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Figure 15. (Color online) Single neutron levels 2p and 1f near the Fermi surface in 
2 


the canonical basis, the rms radii (r)n1;, occupation probabilities v=, j» contributions 
v2 j x (r)nij to the neutron rms radii, and the total neutron rms radii Ry for 
N = 28, 30, 32, and 34 as a function of average pairing gap A. The Fermi surfaces 
for the neutron numbers are plotted as dashed lines. Taken from Ref. [198]. 


In Fig. 15, the single neutron levels ¢; in the canonical basis [188], the rms radii 


(r}nj, the occupation probabilities v2,;, the contributions vf}; x (r}nıj to the neutron 
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rms radius for the orbitals near the Fermi surface (2p and 1f), and the total neutron 
rms radii Ry for four neutron-rich Mg isotopes with N = 28, 30, 32, and 34 are plotted 
as functions of the pairing gap 


1 . 
As N (2j + AW o2. (68) 


nlj 


Here AW) are the diagonal matrix elements of the pairing field in the canonical basis. 

The single neutron levels in canonical basis are obtained by solving the 
corresponding RCHB equations for various pairing strength. In principle, the energy 
levels Eni; = hE? depend on A. However, since this solution was obtained for fixed 
WS potentials, as shown in the upper panels of Fig. 15, these states remain almost 
unchanged with increasing pairing correlations. 

For Mg isotopes with N = 28, 30, 32, and 34, the neutron Fermi surface (shown as 
dashed line in upper panels of Fig. 15) is located just below the two weakly-bound 2p 
levels, between them, crossing the 2pı/2 level, and above the 2p levels. As the pairing 
strength increases, more neutrons are scattered from occupied levels in the Fermi sea to 
empty levels above the Fermi level, and therefore the Fermi surface is raising. 

In Fig. 15, the rms radii of the 2p orbitals are much larger than those of the 
1f orbitals due to the lower centrifugal barrier. With increasing pairing correlations, 
the rms radius decreases for 2p and 1f levels. This is the so-called “pairing anti-halo 
effect” discussed in Ref. [159]. However, this effect concerns only the radii (rm; of the 
individual orbitals. The total neutron radius of the nucleus in Eq. (67) is determined 
by the rms radii of the orbitals and the occupation probabilities, which depend strongly 
on the pairing correlations. 

As shown in Fig. 15, for N = 28 the neutron Fermi surface is just below the two 
weakly-bound 2p levels. Without pairing correlations, the occupation probability is 1.0 
for the 1f7/2 orbital, and it vanishes for the 2p and the 1fs/2 orbitals. As the pairing 
strength increases, the neutrons on 1 f7/2 orbital are scattered to 2p and 1f5/2 orbitals 
which have much larger rms radii. Therefore, the contributions to the total neutron 
rms radius from 2p and 1fs5/2 states grow more than the contribution from 1 f7/2 state 
decreases. As a result, the total neutron rms radius increases monotonically, where the 
contributions from 2p orbitals play a dominant role. 

For N = 30 in Fig. 15, the 2p;/2 and the 1 fs/2 orbitals are empty for zero pairing, 
while the 1 f7/2 orbital is fully occupied and the 2p3/2 orbital is half occupied. With 
increasing pairing strength, the neutrons in the 1f7/2 orbital are scattered to the 2p 
orbitals which contribute strongly to the neutron rms radius. Therefore, the total 
neutron rms radius Ry increases with the growing pairing strength up to A = 0.9 
MeV. Increasing the pairing strength further, the neutrons are continuously to be 
scattered from the 1 f7/2 to the 2p,/2 state, while the neutron number in the 2p3/2 orbital 
decreases slightly and the resonance state 1fs5/2 begins to be occupied. Together with 
the monotonic decreasing rms radii of the individual orbitals, the increase resulting from 
the 2p;/2 and 1fs/2 states is smaller than the decreases of the contributions from the 
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2p3/2 and 1 f7/2 states. Therefore, for very strong pairing correlations, the total neutron 
rms radius finally decreases. As it is clearly seen this decrease is not caused by the 
pairing anti-halo effect, because the radius of the 1f7/2 state stays completely constant. 
It has its origin in the reoccupation caused by pairing. 

For N = 32 in Fig. 15, the 1f7/2 and 2p3/2 orbitals are fully occupied without 
pairing correlations, while the 2p,/2 and 1 fs/2 orbitals are empty. As the pairing strength 
increases, the neutrons in the 1f7/2 and 2p3/2 orbitals are scattered to the 2p1/2 orbital, 
whose rms radius is larger than the ones of 1f7/2 and 2p3/2 levels. The contribution 
of the 2p,/2 orbitals causes an increase in the total neutron rms radius up to a pairing 
gap A = 0.9 MeV. As the pairing strength continues to grow beyond this value, the 
neutrons in the 2p3/2 state begin to be scattered to the 1fs/2 state which provides 
a smaller contribution than the 2p3/2 state, while the occupation probability remains 
almost constant for the 1pı/2 state. Together with the decreasing rms radius of the 
individual orbitals, the total neutron rms radius finally decreases for very large pairing. 

For N = 34, the two weakly-bound 2p levels are fully occupied for zero pairing, 
and the Fermi surface is just above the weakly-bound 2p levels. With increasing pairing 
strength, the neutrons in the 1f7/2 and 2p orbitals are scattered to the 1f5/2 orbital 
with a much smaller rms radius than the 2p orbitals. The contributions to the total 
neutron rms radius from the 1f7/2 and 2p orbitals decrease more than the increase of 
the contribution from 1f5/2 orbital, and therefore the total neutron rms radius finally 
decreases monotonically. Here the two effects, decreasing rms radius of the individual 
orbitals and changes of the occupation probabilities, act in the same direction. The 
decrease of the total neutron radius Ry with increasing pairing correlations is not only 
from the decreasing rms radii of the individual orbitals but also from the change of the 
occupation probabilities. 

It is clear that pairing correlations can change the rms radii of the individual 
weakly-bound orbitals and simultaneously the corresponding occupation probabilities. 
As a result, they have an influence on the nuclear rms radius and the total nuclear 
size. Furthermore, the contributions of the weakly-bound 2p orbitals to the nuclear rms 
radius play an important role. 


5.2. Crucial low l orbitals 


Following the investigation on the impact of nuclear binding on nuclear size, the next 
question is the impact of orbital binding on the nuclear size as a function of the pairing 
potential. In Ref. [198], this problem has been investigated for the orbitals 2p as well- 
bound, weakly-bound, around the continuum threshold, and in the continuum cases. 

In Fig. 16(a), the single neutron levels 1d3/2, 1f7/2, 2p3/2, 2pi/2, and 1fs/2 for 
N = 28 in WS potentials with varying depth such that ¢,,. = —5.0, — 0.5, 0.0, and 
+0.5 MeV are shown. In Fig. 16(b), for each €2p,,,, the change of the total neutron rms 
radius with respect to that of zero pairing, ôRyn = Ry(A) — Ry(A = 0), is plotted as 
a function of the pairing gap. 
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Figure 16. (a) Single neutron levels 1d3/2, 1f7/2; 2p3/2, 2P1/2; and 1fs/2 in WS 
potentials of varying depth with ¢2),,, = —5.0, — 0.5, 0.0, and 0.5 MeV. (b) Changes 
of the total neutron rms radii Ry = Ry(A) — Rn(A = 0) for the various cases of 
E2p,,. aS a function of average pairing gap A. The shadow area represents the results 
calculated with different box sizes in the WS potential for €2p, p= 0.5 MeV. Taken 
from Ref. [198]. 


For €9p,,. = —5.0 MeV which corresponds to a well bound nucleus, the neutron rms 
radius Ry (A) increases slightly with the pairing gap increasing from 0.0 to 1.5 MeV. For 
E2p,/2 = —0.5, 0.0, + 0.5 MeV which correspond to a weakly-bound nucleus, increases 
by more than 0.1 fm are observed. It has to be noticed that for €2p,,, = +0.5 MeV, 
both levels 2p;/2 and 2p3/2 are in the continuum. If the pairing correlation is switched 
on, the system is basically not bound. The eigenvalue problem is solved in a finite box 
and the solution depends on the box size. As an illustration, when the box size changes 
from 20 to 25 fm, the neutron rms radius calculated with a box radius of 25 fm is about 
0.05 fm larger than that obtained with with a radius of 20 fm box for A = 1.5 MeV as 
shown as shadowed region in Fig. 16(b). 

For the discussion on the influence of pairing correlations on the nuclear size so 
far, the mean fields are fixed in the form of spherical WS potentials, and only the 
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pairing strength and the Fermi level for N = 28, 30, 32, 34 are considered. Here, 
with the changing pairing strength and the Fermi level, the single neutron levels stay 
almost constant except slight modification by the pairing field. In Ref. [198], fully self- 
consistent RCHB calculation for neutron-rich nuclei *%424446Mg has been performed, it 
should be emphasized that qualitatively the effect of pairing correlations on the nuclear 
size agrees with that found for fixed WS potentials. 


self-consistent 
calculation 


A (MeV) A (MeV) 


Figure 17. (Color online) Same as Fig. 15, but for “Li and *?Ne. Taken from 
Ref. [198]. 


5.3. Pairing correlations and halos 


In order to investigate the effect of pairing correlations on the development of a nuclear 
halo, self-consistent calculations are performed with different pairing strength for the 
well known neutron halo nuclei Li [34] and **Ne [38], where the low-! orbitals are in the 
continuum [198]. Single neutron levels in the canonical basis, the rms radii (r}mıj, the 
occupation probabilities vz,,, the contributions to the neutron rms radius vz); X (1) nij; 
and the total neutron rms radii are plotted in Fig. 17 as a function of the average pairing 
gap for the nuclei "Li (panel a) and *?Ne (panel b). 
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For "Li in Fig. 17(a), a halo is developed by scattering the neutrons to the 
28/2 orbitals in continuum and the increasing pairing correlations will promote the 
development of this halo. 

For *’Ne in Fig. 17(b), the 2p and 1f;,2 orbitals are in the continuum. Without 
pairing correlations, the last two neutrons occupy the 2p3/2 orbital and therefore this 
nucleus is unbound and the rms radius depends on the box size. Increasing the pairing 
strength, the Fermi surface comes down and becomes negative when the pairing gap 
gets larger than 1.0 MeV. It shows that for the halo nuclei such as 3?Ne, where the low-l 
orbitals are in the continuum and unoccupied without pairing correlation, increasing 
pairing correlations will decrease the size of the nuclear halo, but nonetheless pairing 
plays an essential role: Without pairing the nucleus would stay unbound. 

It is concluded that the weakly-bounded orbitals with low orbital angular momenta 
l in the neighborhood of the threshold to the continuum play an important role [198]. 
The pairing correlations have a two-fold influence on the total nuclear rms radius and 
the nuclear size. First, they can change the rms radii of individual weakly bound 
orbitals and, second, they can change the occupation probabilities of these orbitals. 
With increasing pairing strength, the individual rms radii are reduced (pairing anti- 
halo effect). Meanwhile pairing changes also the occupation probabilities. The total 
nuclear size is determined by the competition between this two effects. 

In Ref. [198], the strength of the pairing correlations was an external variable 
parameter. In realistic nuclei, the strength of pairing depends on the level density 
in the vicinity of the Fermi surface. Therefore, even orbitals with high /-values and 
large centrifugal barriers can contribute indirectly to changes of the nuclear size [38]. If 
they come close to the Fermi surface, they enhance pairing correlations and influence at 
the same time the occupation of low-/ orbitals close to the continuum limit and therefore 
the nuclear size. 


6. Deformed halo and shape decoupling 


Most open shell nuclei are deformed. The interplay between deformation and weak 
binding raises interesting questions, such as whether or not there exist halos in deformed 
nuclei and, if yes, what are their new features. Calculations in a deformed single particle 
model [199] have shown that the valence particles in specific orbitals with low projection 
of the angular momentum on the symmetry axis can give rise to halo structures in the 
limit of weak binding. The deformation of the halo is in this case solely determined by the 
intrinsic structure of the weakly bound orbitals. Indeed, halos in deformed nuclei were 
investigated in several mean field calculations in the past [200, 201, 133]. In Ref. [202], 
it has been concluded that in the neutron orbitals of an axially deformed Woods-Saxon 
potential the lowest-! component becomes dominant at large distances from the origin 
and therefore all 0” = 1/2+ levels do not contribute to deformation for binding energies 
close to zero. Such arguments raise doubt about the existence of deformed halos. In 
addition, a three-body model study [203] suggested that it is unlikely to find halos in 


201707.00914v1 


chinaXiv 


CONTENTS Al 


deformed drip line nuclei because the correlations between the nucleons and those due 
to static or dynamic deformations of the core inhibit the formation of halos. 

Therefore a model which provides an adequate description of halos in deformed 
nuclei must include in a self-consistent way the continuum, deformation effects, large 
spatial distributions, and the coupling among all these features. In addition it should 
be free of adjustable parameters in order to make predictions reliable. 

Over the past years, lots of efforts have been made to develop a deformed 
relativistic Hartree (RH) theory [135] and a deformed relativistic Hartree-Bogoliubov 
theory in continuum (DRHBc theory) [108]. As a first application, halo phenomena in 
deformed nuclei have been investigated within the DRHBc theory [30]. The detailed 
theoretical framework is available in Ref. [31]. The DRHBc theory is useful to 
answer questions concerning exotic nuclear phenomena in nuclei near the drip line 
[200, 199, 204, 140, 203, 201], such as whether there are deformed halos or not and 
what new features can be expected in deformed exotic nuclei. In Ref. [205], the DRHBc 
theory is extended to incorporate the blocking effect due to an odd nucleon. In such 
a way, pairing correlations, continuum, deformation, blocking effects, and extended 
spatial density distributions in exotic odd-A or odd-odd nuclei can be taken into account 
microscopically and self-consistently. In recent years, the CDFT with the density- 
dependent meson-nucleon couplings have attracted more and more attention owing to 
improved descriptions of the equation of state at high density, the asymmetric nuclear 
matter and the isovector properties of nuclei far from stability. Therefore, it’s necessary 
to develop a density-dependent deformed relativistic Hartree-Bogoliubov (DDDRHB) 
theory in continuum. The challenge is the treatment on the density-dependent couplings, 
meson fields, and potentials in axially deformed system with partial wave method. 
Quite recently, the DDDRHB theory in continuum was developed in Ref. [109]. In 
this Section, the halo in deformed exotic nucleus and its mechanism as well as possible 
shape decoupling will be briefly reviewed. 


6.1. Neutron separation energy and nuclear size 


Figure 18 shows the neutron Fermi energy A, and two neutron separation energy 52, of 
Mg isotopes calculated in the DRHBc theory with the parameter sets NL3 and PK1. 
Available data of Sə» taken from [183] are also included for comparison. Except the 
different prediction of the drip line nucleus, the neutron Fermi surfaces and two neutron 
separation energies from both parameter sets are very similar. The calculated two 
neutron separation energies Sn of Mg isotopes agree reasonably well with the available 
experimental values except for °7Mg. The large discrepancy in 3?°Mg is connected to the 
shape and the shell structure at N = 20. For details, see Ref. [31]. 

Experimentally Mg has been observed [206]. Theoretically there are different 
predictions on the last bound nucleus in Mg isotopes, e.g., “Mg from the finite range 
droplet model [207], "Mg from a macroscopic-microscopic model [208], a RMF model 
with the parameter set NLSH [209], and a Skyrme HFB model solved in a 3-dimensional 
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Figure 18. (Color online) The neutron Fermi energy Àn (left panel) and two neutron 
separation energy $2, (right panel) of Mg isotopes calculated by DRHBc theory with 
the parameter sets NL3 and PK1. The data of Sə» (labeled as “Exp.”) are taken 
from [183]. Taken from Ref. [31]. 


Cartesian mesh [210], and Mg from a Skyrme HFB model solved in a transformed 
harmonic oscillator basis [64] and the HFB21 mass table [211]. Therefore the prediction 
of the neutron drip line nucleus in Mg isotopes is model and parametrization dependent. 
In the DRHBc calculations with the parameter set NL3, 4°Mg is the last nucleus of which 
the neutron Fermi surface is negative and the two neutron separation energy is positive. 
However, with the parameter set PK1, Mg is predicted to be the last nucleus within 
the neutron drip line. 

In Fig. 19, the rms radii for Mg isotopes are plotted as functions of the neutron 
number, including the neutron radii Ra, proton radii Rp, matter radii R;,, the ro Al/3 
curve, and the available data for R, [212, 213]. The calculated matter radius follows 
roughly the ro. A'/? curve up to A = 34. From *6Mg on, the matter radius lies much high 
above the roA! curve. This indicates some exotic structure in these nuclei. 


6.2. Density distributions and nuclear shape 


Figure 20 shows neutron density profiles of even-even Mg isotopes with A > 28 
calculated in DRHBc theory with the parameter set PK1 [31]. The spherical component 
of the neutron density distribution is represented by py,,=0(r). The density distributions 
along the symmetry axis z (0 = 0°) and perpendicular to the symmetry axis z (0 = 90°) 
are respectively represented by p(z, rı = 0) and p(z = 0,r_) with ry = yr? 4+ y?. 
With increasing A, the spherical component of the neutron density distribution 
Pn,=0(7) changes quickly at **Mg and abruptly at Mg. The density distribution along 
the symmetry axis pa(z,rı = 0) changes abruptly from **Mg to **Mg. This can be 
understood easily from the shape change in these two nuclei: *?Mg is spherical but *4Mg 
is prolate which means that it is elongated in the z axis. In the direction perpendicular 
to the symmetry axis, the neutron density p,(z = 0,r,) of Mg extends very far away 
from the center of the nucleus and a long tail emerges, revealing the formation of a halo. 
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Figure 19. (Color online) The root mean square (rms) radii for Mg isotopes calculated 
in DRHBc theory plotted as functions of the neutron number. The neutron radii 
Ry, proton radii Rp, matter radii R,, the roA‘/> curve, and the available data for 
R (212, 213] are given. The neutron radius Rn, proton radius Rp and total radius Ry 
of Mg isotopes are calculated with PK1. The ro A!/° curve is included to guide the eye. 
The data for matter radii labeled by “Exp.” are taken from Refs. [212, 213]. Taken 
from Ref. [31]. 


By comparing pn(z,rı = 0) and p,(z = 0,r1) for “Mg, it was found that in the 
tail part, the neutron density distributes more along the direction perpendicular to the 
symmetry axis. Since this nucleus as a whole is prolate, it indicates that the tail part 
has a different shape from the nucleus itself. This is the shape decoupling which was 
first predicted in Ref. [30]. 

The comparison of quadrupole deformation Ø between the deformed RHB 
calculation and the experiment is given in Fig. 21. The experimental value of 8 is 
extracted from the measured B(E£2 : 0f — 27) and only the absolute value could be 
obtained [214]. The results of NL3 are almost the same as the results of PK1 except for 
30Mg. It is slightly oblate with NL3 but it has a prolate shape with PK1. For Mg, the 
gap between the neutron levels 1d3/2 and 1 f7/2 is almost 7 MeV which results in a strong 
closed shell at N = 20. Therefore the deformed RHB calculations with both parameter 
sets predict spherical shapes for this nucleus. This results in the large discrepancy in 
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Figure 20. (Color online) Neutron density profiles of even-even Mg isotopes with 
A > 28 calculated in DRHBc theory with the parameter set PK1. The upper panel 
represents the spherical component of the neutron density distribution pn,\à=0(r). The 
middle panel refers to the density distribution py(z,r,_ = 0) along the symmetry axis 
z (6 = 0°) with ry = yx? + y? and the lower panel refers to pn(z = 0,r1) that 
perpendicular to the symmetry axis z (0 = 90°). Taken from Ref. [31]. 
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Figure 21. (Color online) The quadruple deformation parameter 6 of Mg isotopes 
calculated with the parameter sets NL3 and PK1. The data labeled with “Exp.” are 
taken from [214]. Taken from Ref. [31]. 


So, for **Mg as is seen in Fig. 18. Other mean field models also predict spherical or 
almost spherical shapes for **Mg [215, 209, 210, 216, 217]. However, this nucleus is well 
deformed as is indicated by the large B(E2 : Of — 27) value [218]. The same is true 
for some neighboring nuclei with N ~ 20. These nuclei form the so called “island of 
inversion” [219, 220, 221, 222] which is related to the quenching of the N = 20 shell 
closure. In order to describe the shell quenching effect at N = 20, one has to go beyond 
the mean field [223, 217]. 


6.3. Deformed halo and shape decoupling 


The deformed relativistic Hartree-Bogoliubov theory in continuum with the parameter 
set PK1 predicted the last bound neutron rich nucleus for Mg isotopes as **Mg [31]. 
There are both prolate and oblate stable solutions for *7Mg. The ground state of Mg 
is prolate. With the triaxially deformed relativistic mean field calculations [224], it is 
shown that the obtained oblate minimum is not an isomeric state but a saddle point in 
the potential energy surface. 

It is found that the ground state of Mg is well deformed with 8 ~ 0.41 and two 
neutron separation energy Son © 0.22 MeV [31]. The density distribution of such weakly 
bound nucleus has a very long tail in the direction perpendicular to the symmetry axis, 
see Fig. 20, which indicates the prolate **7Mg has an oblate halo. 

The density distribution is divided into contributions of the oblate “halo” and of 
the prolate “core” in Fig. 22. The density distribution of this weakly bound nucleus has 
a very long tail in the direction perpendicular to the symmetry axis which indicates the 
prolate nucleus *?Mg has an oblate halo and there is a decoupling between the shapes 
of the core and the halo. It should be noted that in the calculations with NL3, “Mg is 
the last nucleus [30]. 

The single particle spectrum around the Fermi level for the ground state of Mg 
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Figure 22. (Color online) Density distributions of the ground state of 4?7Mg with the 
z axis as the symmetry axis: (a) the neutron halo, and (b) the neutron core. Taken 
from Ref. [31]. 


is shown in Fig. 23 [31]. The good quantum numbers of each single particle state are 
also shown. The occupation probabilities v? in the canonical basis have BCS-form [188] 
and are given by the length of the horizontal lines in Fig. 23. The levels close to the 
threshold are labeled by the number 7 according to their energies, and their conserved 
quantum number Q” as well as the main spherical components are given at the right 
hand side. The neutron Fermi level is within the pf shell and most of the single particle 
levels have negative parities. Since the chemical potential An is close to the continuum, 
orbitals above the threshold have noticeable occupations due to the pairing correlations. 
The single neutron levels of Mg can be divided into two parts, the deeply bound levels 
(Ecan < —2 MeV) corresponding to the “core”, and the remaining weakly bound levels 
close to the threshold (€¢a, > —0.3 MeV) and in the continuum corresponding to the 
“halo”. 

As discussed in Refs. [30, 31, 225], the shape of the halo originates from the intrinsic 
structure of the weakly bound or continuum orbitals. By examining the neutron density 
distribution, it was found that for the ground state of “?Mg, the halo is mainly formed 
by level 4 and level 5. A detailed analysis of the intrinsic properties revealed that in 
42Mg the shape of the halo is oblate and decouples from the prolate core. 

We note that this kind of shape decoupling was later also found in non-relativistic 
HFB calculations [226]. An exotic “egg” -like structure consisting of a spherical core plus 
a prolate halo was predicted in **Ne, in which the near-threshold nonresonant continuum 
plays an essential role. 


6.4. Recent progress 


For odd particle system, the halo could depend strongly on the interplay among the 
odd-even effects, continuum and pairing effects, deformation effects, etc. Combining 
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Figure 23. (Color online) Single neutron levels of ground state of 4*Mg in the 
canonical basis as a function of the occupation probability v?. The blue dashed line 
corresponds to the BCS-formula with an average pairing gap. Taken from Ref. [31]. 


the necessity for the improved descriptions of the equation of state at high density, 
the asymmetric nuclear matter and the isovector properties of nuclei far from stability 
with the density-dependent meson-nucleon couplings [227, 228, 229, 230, 138, 182], it 
is necessary to extend the DRHBc theory with the density-dependent meson-nucleon 
couplings. Here we will briefly review these progress, i.e., the exotic nuclear structure in 
unstable odd-A or odd-odd nuclei [205] and the density-dependent deformed relativistic 
Hartree-Bogoliubov (DDDRHB) theory in continuum [109]. 

In order to describe the exotic nuclear structure in unstable odd-A or odd-odd 
nuclei, the blocking effect of one or several nucleons has to be taken into account. In 
Ref. [205], the DRHBc theory [108, 30, 31] is extended to incorporate the blocking effect 
due to odd particle system. In such a way, pairing correlations, continuum, deformation, 
blocking effects, and extended spatial density distributions in exotic odd-A or odd-odd 
nuclei can be taken into account microscopically and self-consistently. 

For an odd nucleon system with the i-th level blocked, the time reversal symmetry 
is violated and there appear currents in the system. In principle, one has to diagonalize 
the RHB equation twice. Using the equal filling approximation which is usually 
made [231, 232], the currents are neglected and the two configurations of a particle in the 
time-reversal partner space are averaged in a statistical manner [205]. The corresponding 
currents cancel each other and in this way, in each step of the iteration meson fields 
with the time reversal symmetry is obtained. In practice, the density matrix p and the 
abnormal density « in two subspaces are averaged and replaced by 


1 x x 
p = Puxm ag 2 (UnU _ Vie Va) , (69) 
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Figure 24. (Color online) Neutron single particle levels in the canonical basis and 
their occupation probability v? for 12O from the DRHBc and RCHB calculations. 
Taken from Ref. [205]. 


K = Kuxm — : (UVa T Viz UE.) ? (70) 
where Vp, and Up, are column vectors in the matrices V and U corresponding to the 
blocked level. 

Using the density functional PK1 [138] and a zero range density dependent pairing 
force [39, 31], the deformed RHB equations are solved in a spherical Dirac Woods- 
Saxon basis [65] and the deformed potentials and densities are expanded in terms of the 
Legendre polynomials P\(cos@) (for details see Refs. [108, 30, 31]). In order to check 
the accuracy of the DRHBc code with the blocking effect, the spherical nucleus 19O is 
investigated and the results of the DRHBc code allowing only the spherical components 
of the fields, i.e., A = 0, are compared with those obtained in RCHB theory [39]. It is 
found that both calculations agree well with each other with high precision [205]. 

In Fig. 24, the neutron single particle levels in the canonical basis for 19O is shown 
in comparison with the RCHB results in order to understand the details and validity 
of the DRHBc calculation. The length of each level is proportional to the occupation 
probability. In the RCHB calculation, the blocked orbital is 1d5/2. In the DRHBc code, 
a spherical solution is enforced by allowing only 4 = 0 component in the Legendre 
expansion mentioned earlier and therefore the three sublevels of the 1d5/2 orbital with 
Q7 = 1/2t, 3/27, and 5/2+ are degenerate. The present results in Fig. 24 are obtained 
with the level Q7 = 1/2+ blocked in the DRHBc code, but blocking any of these three 
levels gives the same results. From Fig. 24 one finds nice agreement between the results 
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Figure 25. (Color online) (a) Neutron density profiles of ‘tC and 15C in DRHBc 
calculations with the parameter set PK1. (b) Neutron single particle levels in the 
canonical basis and their occupation probability v? for 15C in DRHBc calculation. 
Taken from Ref. [205]. 


from the DRHBc and RCHB calculations. 

An application of the DRHBc theory for an exotic and deformed nucleus with an 
odd number of particles is given in Ref. [205] for the deformed halo candidate nucleus 
15C [233] with the parameter set PK1 [138]. In contrary to the spherical even-even 
nucleus 14C, the DRHBc theory predicts a deformation 3 = 0.25 for the odd-A nucleus 
15C, The neutron root mean square radii of C and °C are 2.56 fm and 2.79 fm, 
respectively. 

The DRHBc calculations for °C lead to a deformed density distribution with an 
axial symmetry. In Fig. 25, the neutron densities along the symmetry axis p,(rz,r = 0) 
and perpendicular to the symmetry axis p,,(r, = 0,r_) are plotted as dashed and dashed- 
dotted lines. The spherical density for 14C is also included as a reference by a solid line. 
It is interesting to see that, in the direction perpendicular to the symmetry axis, the 
neutron densities are almost the same for 14C and °C at least when rı < 6 fm. Along 
the symmetry axis the neutron density of °C extends much further than that of 14C. 
This is partly due to that !°C is prolate and the weakly bound 1/27 level is occupied. 

The single particle levels of °C in the canonical basis are plotted in Fig. 25. Since 
it is a deformed nucleus, the 1p3/2 orbital is split into two levels with Q7 = 1/27 and 
Q7 = 3/27, respectively. There is one neutron occupying in Q7 = 1/2% level near the 
threshold. The occupation probability v? = 0.5, indicating that it is averaged over the 
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two configurations with Q7 = +1/2*. Because of the deformation, this level is a mixture 
of the spherical orbitals 1ds/2 (62%) and 281). (36%). The weakly-bound feature and 
the relatively large s-wave component of this level results in that the neutron density of 
15C extends further along the symmetry axis. 

After the comparison between the DRHBc and the spherical RCHB results for °O 
as well as the application for the deformed neutron-rich nucleus °C, the DRHBc with 
blocking effect proved to be a powerful tool for the study of the odd-A or odd-odd exotic 
nuclei. 

In the development of the DRHBc with the density-dependent meson-nucleon 
couplings [109], the challenge is the treatment of the density-dependent couplings, meson 
fields, and potentials in axially deformed system with partial wave method. 

In Fig. 26, by taking the o meson as an example, the partial waves of 
density-dependent coupling strengths g,,,(r) for **Mg are shown with golpe) = 
Xa 9c,,(")P\(cos 6). It can be seen that the major component is that with A = 0, 
which is more than one order of magnitude larger than the others. The amplitudes of 
Jo,. decrease quickly with increasing À and become negligible when À > 8. This means 
that a cutoff Ag max = 10 is enough to provide accurate results. 


r (fm) 


Figure 26. (Color online) Decompositions of coupling strength g,, with A = 0 (a) 
and \ = 2,4,6,8,10 (b) for °°Mg calculated with density-dependent DRHBc theory. 
Taken from Ref. [109]. 


201707.00914v1 


chinaXiv 


CONTENTS 51 


By comparing the calculated bulk and single particle properties of neutron-rich 
nucleus **Mg in spherical case with Ap,max = 0 with those obtained by RCHB [39], 
the density-dependent DRHBc theory is verified [109]. It therefore opens the door to 
investigate the deformed exotic nuclei with effective density-dependent interactions. 


7. Summary and perspectives 


In this Topical Review, we have reviewed the progresses on the covariant density 
functional theory in continuum for the neutron halo phenomenon in both spherical 
and deformed nuclei. 

In the mean field study of nuclear halos, a proper description of nuclear pairing and 
the continuum contribution is very important. To include properly the contribution of 
continuum states, the Bogoliubov transformation has been justified to be very useful. 
The relativistic continuum Hartree-Bogoliubov (RCHB) theory has been developed in 
which the coupled differential equations are solved in coordinate (r) space and the mean 
field effects of the coupling to the continuum can be fully taken into account. 

From the RCHB theory, a self-consistent description of the neutron halo in 'Li 
has been achieved. The calculated density distribution agrees well with the experiment. 
It has further been shown that the halo in "Li is formed by Cooper-pairs scattered 
in the two levels 1pı;2 and 2s,/2 with the former below the Fermi surface and the 
latter in the continuum. In the RCHB theory, a halo appears if there are several single 
particle valence orbitals with small orbital angular momenta and correspondingly small 
centrifugal barriers close to the continuum limit. 

For medium and heavy nuclei, the RCHB theory has treated the influence of 
pairing correlations and many-body effects properly and predicted the existence of more 
neutrons in the halo. Giant neutron halos, consisting of up to six halo neutrons, have 
been predicted in Zr and Ca isotopes close to the neutron drip line. When approaching 
the neutron drip line, the two-neutron separation energies for several nuclei in Ca and 
Zr isotopes are extremely close to zero, which means that the valence neutrons gradually 
occupy the loosely bound states and the continuum. 

The phenomena of halo and giant halo of neutron-rich even-Ca isotopes were 
investigated and compared in the framework of the RCHB and non-relativistic Skyrme 
HFB calculations. It has been found that although halo phenomena exist for Ca isotopes 
near neutron drip line in both calculations, the halo in the Skyrme HFB calculations 
starts at a more neutron-rich nucleus than that in the RCHB calculations, and the RCHB 
calculations have larger neutron root-mean-square radii systematically in N > 40 than 
the Skyrme HFB calculations. The former difference comes from the difference in shell 
structure. The reasons for the latter can be partly explained by the neutron 3s1/2 orbital, 
which causes more than 50% of the difference for neutron radii at Ca. 

The relativistic Hartree-Fock-Bogoliubov theory in continuum (RHFBc theory) 
with density-dependent meson-nucleon couplings has been developed. With this theory, 
the influence of exchange terms on the formation of neutron halos has been studied. 
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Giant halos as well as ordinary ones were found in Ce isotopes close to the neutron drip 
line. It was demostrated that the shell structure evolution, when approaching the drip 
line, is important for the appearance of neutron halos. 

The RCHB theory has also been extended to including hyperons and even Ca 
isotopes ranging from the proton drip line to the neutron drip line have been studied 
systematically. Compared with ordinary nuclei, the additional one or two A-hyperons 
lower the Fermi level and the last bound hypernuclei have two more neutrons than the 
corresponding ordinary nuclei. Giant halo phenomena due to the pairing correlation and 
the contribution from the continuum have been predicted in Ca hypernuclei, similar to 
those appearing in ordinary Ca isotopes. 

A deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc theory) 
has been developed aiming at a proper description of exotic nuclei, particularly those 
with a large spatial extension. In order to give an adequate consideration of both the 
contribution of the continuum and the large spatial distribution in exotic nuclei, the 
deformed RHB equations were solved in a Woods-Saxon (WS) basis in which the radial 
wave functions have a proper asymptotic behavior at large distance from the nuclear 
center. 

Halo phenomena in deformed nuclei were investigated within the DRHBc theory. 
These weakly bound quantum systems present interesting examples for the study of 
the interdependence between the deformation of the core and the particles in the halo. 
The Mg and Ne isotopes have been studied in details. It was found that in 4?*4Mg, 
there appears a decoupling of the halo orbitals from the deformation of the core: The 
core is prolate, but the halo has a slightly oblate shape. The generic conditions for the 
occurrence of this decoupling effects have been given: the existence and the deformation 
of a possible neutron halo depends essentially on the quantum numbers of the main 
components of the single particle orbitals in the vicinity of the Fermi surface. 

Definitely there exist lots of open questions on the study of nuclear halos with the 
covariant density functional theory. We list several of them below. 

In the RCHB, RHFBc and DRHBc theories, the box boundary condition has been 
applied and the continuum is discretized. Therefore one cannot get the width for 
resonances. With the Green’s function method, the correct asymptotic behavior can be 
given on the wave function for the continuum states. The Green’s function method has 
been implemented in the non-relativistic HFB theory [55]. The effort to implement the 
Green’s function approach in the RHB theory is undertaken and promising preliminary 
results have been obtained [58]. 

In the systematic study of giant halos in severl isotope chains, it has been found 
that the appearance of giant halos and the contribution of continuum are very important 
in determining the boundary of nuclear chart. Furthermore, the deformation effects will 
also influence very much the position of drip lines. More systematic study is now carried 
out to reveal the global impact of the continuum and deformation effects on the nuclear 
boundary. 

If the core and the halo of a nucleus have different shapes, many interesting 
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questions arise: For example, What will happen when the nucleus rotates? What 
role could the shape fluctuation play? To answer these questions, future works are 
welcome to implement the angular momentum projection in the DRHBc theory and to 
use the generator coordinate method (GCM) or to solve Bohr Hamiltonian based on the 
potential energy surface obtained from constraint DRHBc calculations. 

The possible experimental signals of shape decoupling in deformed halo nuclei 
are certainly interesting topics. The halo feature is connected with relatively large 
cross sections and narrow longitudinal momentum distributions in knockout reactions. 
The decoupling between the deformations of the core and the halo may manifest itself 
by some new experimental observables, e.g., the double-hump shape of longitudinal 
momentum distribution in single particle removal reactions and new dipole modes, etc. 
In particular, a combination of the experimental method proposed in Ref. [234] and the 
theoretical approach developed in Ref. [235] would be useful in the study of longitudinal 
momentum distribution in single particle removal reactions with deformed halo nuclei as 
projectiles. The shape decoupling effects may also has some influence on the sub-barrier 
capture process in heavy ion collisions [236]. 
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